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Abstract 

In this paper, we extend the uniform regularity estimates obtained by M. Avellaneda and F. Lin in GE 
to the more general second order elliptic systems in divergence form {C e ,E > 0}, with rapidly oscillating 
periodic coefficients. We establish not only sharp W 1,p estimates, Holder estimates, Lipschitz estimates 
and non-tangential maximal function estimates for the Dirichlet problem on a bounded C 1,v domain, but 
also a sharp O(e) convergence rate in Hq (12) by virtue of the Dirichlet correctors. Moreover, we define 
the Green’s matrix associated with C e and obtain its decay estimates. We remark that the well known 
compactness methods are not employed here, instead we construct the transformations (11.111) to make 
full use of the results in m- 


1 Introduction and main results 

The main purpose of this paper is to study the uniform regularity estimates for second order elliptic 
systems with lower order terms, arising in homogenization theory. More precisely, we consider 

C e = —div [A{x/e)V + V {x/e)\ + R(x/e)V + c(x/e) + XI, 

where A is a constant, and I = ( S a denotes the identity matrix. In a special case, let A = I = 1, V = B, 
c = 0, and W = div(U), the operator C e becomes 

£ e = -A + ^W(x/e) + A, 

where W is the rapidly oscillating potential term (see [7, pp.91]). It is not hard to see that the uniform 
regularity estimates obtained in this paper are not trivial generalizations of mm, and they are new even for 
£ e - 

Let 1 < i,j < d, 1 < a, /3 < m, where d > 3 denotes the dimension, and m > 1 is the number of equations 
in the system. Suppose that the measurable functions A = (aff) ■ -A M m2x<i2 , V = ( V'^) : -A M m2xd , 

B = (B^) : R. d -A M m2xd , c = (c a/3 ) : -A M™ 2 satisfy the following conditions: 

• the uniform ellipticity condition 

m 2 < for y € and e = (£f) G where /x > 0; (1.1) 

(The summation convention for repeated indices is used throughout.) 
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• the periodicity condition 

A(y + z) = A(y), V(y + z) = V(y), B(y + z) = B(y), c{y + z) = c(y) for y G R d and z G Z d ; (1.2) 

• the boundedness condition 

m ax{||V r || i oo (R d ) , ||c|| L oo (Rd) } < m , where ki>0; (1.3) 

• the regularity condition 


Set k = max{/ci, /C 2 }, and we say A G A(/r,r, k ) if A = A{y) satisfies conditions (II. ID . (11.21) and CUD- 
Throughout this paper, we always assume that fl is a bounded C 1,r> domain with 77 G [r, 1), and L e = 
—div[A(x/e)V] is the elliptic operator from [3l|6j, unless otherwise stated. 

The main idea of this paper is to find the transformations (11.111) between two solutions corresponding 
to C £ and L e such that the regularity results of L e can be applied to C e directly. Particularly, to handle the 
boundary Lipschitz estimates, we define the Dirichlet correctors <h £i fc = (<f>“^), 0 < k < d, associated with 


C £ as follows: 


for k = 0 , and 


L £ (® £ ,k) = div(W) in LI, 


= I on <90 


£ £ (< fc )= 0 infi, 


®ek = P k on ^Lt 


(1.5) 


( 1 . 6 ) 


for 1 < k < d, where V £ (x) = V{x/e), = (4?^, ■ ■ ■ , ) G and Pf! = Xk(0, • • • , 1, ■ ■ ■ ,0) 

with 1 in the (3 th position. We remark that (11.61) was studied in |3lf38]. but () 1.5 1) has not yet been developed. 
Here we show that <f> £) o ought to be of the form in (11.51) . and its properties are shown in section 4. 

For Neumann boundary conditions, a significant development was made by C.E. Kenig, F. Lin and Z. 
Shen [3T|, where they constructed Neumann correctors to verify the Lipschitz estimates for L e . Recently, 
S.N. Armstrong and Z. Shen [2] found a new way to obtain the same results even without Dirichlet correctors 
or Neumann correctors in the almost periodic setting. We plan to study uniform regularity estimates for C £ 
with Neumann boundary conditions in a forthcoming paper. 

The main results are as follows. 

Theorem 1.1 {W 1,p estimates). Suppose that A G VMO(M d ) satisfies (11.11) . (11.21) . and other coefficients of 
JC £ satisfy £T3D- Let 1 < p < 00 , f = (/“) € L p (Ll] R md ), F G L q (Ll-, R m ) and g € B l ~P p (dLl-, M m ), where 
q = pTd */ P> d^T > an d 9 > 1 */1 < p < d~[ ■ Then the Dirichlet problem 


C £ (u £ ) = div(/) + F in LI, 
u £ = g on dLl 


(1.7) 


has a unique weak solution u £ € VF 1 ’ p (H;M m ), whenever A > Ao and Ao = Ao(/r, k, m, d) is sufficiently large. 
Furthermore, the solution satisfies the uniform estimate 


||V« £ ||iP(n) < C{\\f\\ L p(n) + ||T|| Li{ f2) + 
where C depends on p,ui(t), k, X,p,q,d,m and LI. 

Note that A G VMO(M d ) if A satisfies 

f A(y) — + A dy < u>(t), and lim u>(t) = 0, 
Jb(x,p) Jb(x.p) *-»0+ 


( 1 . 8 ) 


sup 

xeRd JB(x,p) 

0 <p<t 







3 


and B" ,p (<912; M m ) denotes the L p Besov space of order a (see [T]). We mention that for ease of notations 
we say the constant C depends on ui instead of u:(t) in the rest of the paper. We will prove Theorem II.II in 
Section 3 by using bootstrap and duality arguments. We mention that there are no periodicity or regularity 
assumptions on the coefficients of the lower order terms in Theorems 11.11 and 11.21 The estimate ([1.8D still 
holds when 12 is a bounded C 1 domain (see TO- 

Results of the W 1,p estimates for elliptic or parabolic equations with VMO coefficients can be found 
in |101tl2iri61l34lf35| . In the periodic setting, similar estimates for parabolic systems, elasticity systems, and 
Stokes systems were obtained by mmm, respectively. Also, the uniform W l,p estimates for L e with 
almost periodic coefficients were shown in [2] recently. 

Theorem 1.2 (Holder estimates). Suppose that the coefficients of C £ satisfy the same conditions as in 
Theorem fLll Let p > d, f = (/“) € L p (12;M md ), F G L q (Q; M m ), and g G C°’ CT (<912; M m ), where q = 
and <7 = 1 — d/p. Then the weak solution u £ to (11.71) satisfies the uniform estimate 

ll^ellc'O.'TQ) < C{\\f\\ L P(Cl) + ll-F||z,®(n) + IMIc 0 ^ (dCl )} > (1-9) 


where C depends on g, w, k, A, p, a, d, m and 12. 

The estimate (11.91) is sharp in terms of the Holder exponent of g. If g G C' 0,1 (<912; M m ), (11.91) is just 
Corollary 13.81 The uniform Holder estimates for L e were given in [3] by the compactness method which also 
works for non-divergence form elliptic equations (see a). However, we can not derive the sharp estimate 
by simply applying this method. So we turn to study the Green’s matrix G e {x,y) associated with C £ and 
obtain the decay estimates 


\Qe(x,y)\ < 


C 

x — y\ d ~ 2 


min 



d% 

x — y\ a 


d a ' d a d a ' 

Liy LL X Uy 

x-y/ 7 " \x-y\ a + a ' y 


V x,y £ 12, x^y, 


where a. a’ G (0,1), d x = dist(x, 512) denotes the distance between x and <912, and C is independent of e (see 
Theorem 13.lip . Then we prove Theorem 11.21 through a subtle argument developed by Z. Shen [40], where 
he proved a similar result for L e in the almost periodic setting. 

The existence and some related properties of the Green’s matrix with respect to L\ were studied by S. 
Hofmann and S. Kim [26]. We also refer the reader to m for parabolic systems, and [23][36] for the scalar 
case. 


Theorem 1.3 (Lipschitz estimates). Suppose that A G A (h,t,k), V satisfies (jl.2l) . (11.41) . B and c satisfy 
(11.31) . and A > Ao- Let p > d and 0 < a < rj. Then for any f G C 0,<T (12;M mrf ), F G L p (12;M m ), and 
g G C' 1,fT (512; M m ) ; the weak solution to (|1.7I) satisfies the uniform estimate 

I|V« £ ||loo(q) < C {||/||c°' CT (n) + H-^llLP(f 2 ) + llfl'llc 1 ’ CT ( 0 f 2 )} > (1-10) 

where C depends on g, r, k, X,p, d, m, a, g and 12. 

The estimate (I1.10|) can not be improved even if the coefficients of C £ and 12 are smooth, since the 
corrector xo defined in ()2.1I) is a counter example. Here we use two important transformations 

u £ = [I + £Xo(x/e)\v £ , and u £ = § £ $v £ (1.11) 


to deal with the interior and global Lipschitz estimates, respectively. We explain the main idea as follows: 


(Di) 


C £ {u £ ) = div(/) + F 
u £ = g 


in 12, 
on <912. 


Ue = $£, 0 ^ 
- > 


(D 2 ) 


L £ {v £ ) = div(/) + F 
Ve = 9 


in 12, 
on <912. 


Note that <h £i o is not periodic, which is the main difficulty to overcome. So we rewrite (Di) as (D 2 ) to keep 
L £ periodic, while the price to pay is that the new source term / involves Vu £ . As we mentioned before, 
there is no uniformly bounded Holder estimate for Vu £ . Fortunately, it follows from Theorem 11.21 that 


V < & £ ,o||c'°’ ir 1 (n) — 0( £ CT1 ) an d 11V i/gr 11 c°,o-i(j^) — 0(s a2 ) as e —» 0, 


( 1 . 12 ) 
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where 0 < o\ < 02 < 1 are independent of e (see Lemma 14.91 and 14.101) . Together with an important 
consequence of Lemma 14.81 

|| < f > e,o ~~ I\\L°°{n) = 0 {e) as e — > 0, (1-13) 

we obtain that / is uniformly Holder continuous through the observation that the convergence rate in (| 1. 13 jl 
is faster than the divergence rate in (11.121) as e —>■ 0. Also, Theorem 11.11 implies F £ L p [Ll) with p > d. 
Thus we can employ the results in [3] immediately, and the proof of Theorem 11,31 is finalized by a suitable 
extension technique. 

We remark that the compactness argument for our elliptic systems is also valid, however it would be 
much more complicated. For more references, C. E. Kenig and C. Prange [32] established uniform Lipschtiz 
estimates with more general source terms in the oscillating boundaries setting, and the same type of results 
for parabolic systems and Stokes systems were shown in [dEI], respectively. 

Theorem 1.4 (Nontangential maximal function estimates). Suppose that A £ A(//, r, k), V,B satisfy (11.21) 
and m, c satisfies dm and A > Ao- Let 1 < p < 00 , and u e be the solution of the L p Dirichlet 
problem C £ {u £ ) = 0 in LI and u £ = g on dLl with (u e )* £ L p {dLl), where g € L p {dLl]W n ) and (u e )* is the 
nontangential maximal function. Then 


IK«e)*lli>(an) < CplMlLP(an)) (1-14) 

where C p depends on y,r, k, X, d,m,p,p and LI. Furthermore, if g £ L 00 (dLl]M rn ), we have 

11 111/°° (fi) < C1bllz,°°(afi)j (1-15) 

where C depends on y,r, k, A, d,m,g and H. 

The estimate (11.151) is known as the Agmon-Miranda maximum principle, and (u e )* is defined in (14.401) . 
We remark that the proof of Theorem 11.41 is motivated by [31EI1S8] • Define the Poisson kernel associated 
with C £ as 

Vf(x,y) = -n j (x)aff(x/e)-^{G^(x,y)} -n j {x)B^{xfe)gf 1 {x,y), 

where nj denotes the j th component of the outward unit normal vector of dfl. Due to Theorem 11.31 we 
obtain lV x V y Q £ (x, y)\ < C\x — y\~ d for x,y £ D, and x y (see Lemma 11 .11 1) . which implies the decay 
estimate of V £ (see ([4.371) 1. Thus the solution u £ can be formulated by ([4.381) . Note that V £ is actually 
closely related to the adjoint operator C* (see Remark 12.31) . That is the reason why we additionally assume 
m and (11.41) for B in this theorem. We refer the reader to Remark 14.121 for more references on Theorem 

Ol 

Theorem 1.5 (Convergence rates). Let Lt be a bounded C 1,1 domain. Suppose that the coefficients of C £ 
satisfy the same conditions as in Theorem 11.31 and B,c additionally satisfy the periodicity condition (|1.2j) . 
Let u £ be the weak solution to C £ {u E ) = F in LI and u £ = 0 on dLl, where F £ L 2 (D;M m ). Then we have 

a„.f 3 

\\u £ -$ £i0 u-{^ k -pP)— || H i (n) <Ce\\F\\ L 2 [n) , (1.16) 

where u satisfies Cq{u) = F in Ll and u = 0 on dLl. Moreover, if the coefficients of C £ satisfy (ED - ED> 
then 

\\u £ - u\\li(q.) < Ce\\F\\ LP (£i) (1-17) 

holds for any F £ L P (D; M m ), where q = if 1 < p < d, q = 00 ifp > d, and C depends on y, r, k, A, m , d,p 
and Ll. 
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We mention that the estimates (11.161) and (11.171) are sharp in terms of the order of e. The ideas in 
the proof are mainly inspired by [29(30]. It is easy to see \\u £ — u ||x, 2 (fi) = 0(e) is a direct corollary of 
(11.161) or (11.171) . In the case of p = d, (11.171) is shown in Remark 15.31 Moreover, in the sense of “operator 
error estimates” the convergence rate like (11.171) can also be expressed by WCj 1 — C^Wlp^li < Ce, where 
|| ■ is referred to as the ( L p —> L 9 )-operator norm. 

The convergence rates are active topics in homogenization theory. Decades ago, the L 2 convergence 
rates were obtained in [71128] for scalar cases due to the maximum principle. At the beginning of 2000’s, the 
operator-theoretic (spectral) approach was successfully introduced by M. Sh. Birman and T. A. Suslina [5115] 
to investigate the convergence rates (operator error estimates) for the problems in the whole space M rf . They 
obtained the sharp convergence rates 0(e) in the ( L 2 —> T 2 )-operator norm and (L 2 —> H 1 )- operator norm 
for a wide class of matrix strongly elliptic second order self-adjoint operators, respectively. These results were 
extended to second order strongly elliptic systems including lower order terms in [44] . Recently, C. E. Kenig, 
F. Lin and Z. Shen [29] developed the L 2 convergence rates for elliptic systems on Lipschitz domains with 
either Dirichlet or Neumann boundary data by additionally assuming regularity and symmetry conditions, 
while T.A. Suslina [421143] also obtained similar results on a bounded C' 1,1 domain without any regularity 
assumption on the coefficients. We refer the reader to miiai20l[37l[46| and references therein for more 
results. 

In the end, we comment that the above five theorems are still true for d= 1,2. Since we usually have a 
different method to treat the cases d> 3 and d = 1,2 (for example, see [26, pp.2]), we omit the discussion 
about the cases of d = 1,2 here. 

2 Preliminaries 

Define the correctors Xk = (X}f)i 0 < k < d, associated with C £ as follows: 

L\(Xk) = div(V) in M rf , 

Xk G Hl er (Y;R m2 ) and Xkdy = 0 ^ ) 

for k = 0, and 

' iifef + H) = » 

’ xfe-gUr; M m ) 

for 1 < k < d, where Y = [0, l) d = M d /Z d , and Hp er (Y]W n ) denotes the closure of C^ r (Y;W n ) in 
H 1 (Y;W n ). Note that C^^Y^M 771 ) is the subset of C°°(y;M m ), which collects all y-periodic vector-valued 
functions (see m pp.56]). By asymptotic expansion arguments, we obtain the homogenized operator 

Co = -div(AV + V) + BV + c + XI, 
where A = (a°j 3 ), V = (V) a/3 ), B = (B^) and c = (c a/3 ) are given by 



The proof is left to readers (see 0 pp.103] or (28] pp.31]). 

Remark 2.1. It follows from the conditions (11.11) and (11.21) that ^t|£| 2 < < /4|£| 2 holds for any 

^ = (4°) e K md , where /xi depends only on /x. Moreover, if a^ 3 = a^“, then /x|^| 2 < ^ < /U _1 |^| 2 

(see 0 pp.23]). This illustrates that the operator Co is still elliptic. 


m 


and / Xk d V = 0 


( 2 . 2 ) 
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Remark 2.2. We introduce the following notations for simplicity. We write Xk.fix) = Xk{x/e ), A e (x ) = 
A(x/e), V e {x) = V{x/e), B e (x ) = B(x/e), c £ (x) = c(x/s), and their components follow the same abbreviated 
way. Note that the abbreviations are not applied to or ^> £ ^{x). 

Let B = B(x,r ) = B r (x), and KB = B(x,Kr) denote the concentric balls as K > 0 varies, where 
r < 1 in general. We say that Ll is a bounded C 1,T1 domain, if there exist ro > 0, Mq > 0 and {Pi : i = 
1,2, • • • , no} C dLl such that dLl C U^ 1 R(Pj, ro) and for each i, there exists a function G C 1 ,? 7 (]R d_1 ) and 
a coordinate system, such that !?(/}, Coro) fl fi = P(0, Coro) fl {{x',Xd) G : x' G and Xd > fifix')}, 

where Cq = 10(Mq + 1) and fii satisfies 


V’i(O) = 0, and ||V’i|lc 1 >’7(M d - 1 ) < M o- (2.4) 

We set D(r) = D{r,fi) = {(xfxf) G : |x'| < r and ip(x') < Xd < x ') + Cor} and A (r) = A(r, fi) = 
{(a/, T/>(a/)) G : |a/| < r}. In the paper, we say the constant C depends on If, which means C involves 
both Mo and |O|. This is especially important when we do near boundary regularity estimates. Here |S2| 
denotes the volume of LI. We also mention that for any £ C SI, we write Je = f{x)dx = f E f(x)dx , 
and the subscript of f e is usually omitted. 


Remark 2.3. Let C* £ be the adjoint of C £ , given by 


+ B f(- X / £ ) V e} + + ca ^ x / £ K + Auf 


dv\ 


dxi ' Bj 

Then we define the bilinear form associated with C £ as 

J 


dv £ 


B«K,0] = j a {«“i^ + Vffnlf-fi-dx + f a | + cf4 + Xu«}rdx, 


yapdVe_ a/3 a 

dXi ^ £ e + £ 


0 


and the conjugate bilinear form with respect to C* £ as 

b-[,. m- f 1 n'vC* , 

el 01 “ Ja v“«’« + 1 ax , dx + Jn 

for any u £ ,v £ ,fi G ILo(H;E m ). It follows that B e [u e ,u e ] = B *[v £ ,u £ ] and 


< T £ (u £ ),v £ O— 


(ALVu £ + V £ u £ )Vv £ dx - [ u £ div(B £ v £ )dx+ [ (c £ + X)u £ v £ dx 


(2.5) 


=< u £ ,C* e (y £ ) > . 


( 2 . 6 ) 


Lemma 2.4. Let Ll be a Lipschitz domain. Suppose that A satisfies the ellipticity condition m, and other 
coefficients of C £ satisfy (11.31) . Then we have the following properties: for any u,v G f7o(fi;M m ), 

|B £ [u,u]| < C||u|| H i (a) ||u|| H i (a) , and c 0 |M|^i (n) < B e [u, u), whenever A > A 0 , (2.7) 


where Ao = Ao (n,K,m,d) is sufficiently large. Note that C depends on pL, k, X,m,d,Ll, while Co depends on 
pL, k, m, d, Ll. 

Theorem 2.5. The coefficients of C £ and Ao are as in Lemma 12.41 Suppose F G H~ 1 (Ll-,W n ) and g G 
i 72 (< 9 f 2 ;M m ). Then the Dirichlet boundary value problem C £ (u £ ) = F in Ll and u £ = g on dLl has a unique 
weak solution u £ G iL 1 (H), whenever A > Aq, and the solution satisfies the uniform estimate 


ll u e||.ff 1 (n) < C{||F|| H -i (n) + Ibllif 1 / 2 ^)}) ( 2 - 8 ) 

where C depends only on /j,,K,m,d and Ll. Moreover, with one more the periodicity condition on the 
coefficients of C £ , we then have u £ —*• u weakly in i7 1 (fl;M m ) and strongly in L 2 (H;M m ) as e — > 0, where u 
is the weak solution to the homogenized problem Cfiu) = F in LI and u = g on dLl. 
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Remark 2.6. The proof of Lemma 12.41 follows from the same argument in the scalar case (see |15]l'22j). 
Theorem 12.51 involves the uniqueness and existence of the weak solution to (11.71) , and the so-called homog¬ 
enization theorem associated with C £ . The proof of Theorem 12.51 follows from Lemma 12.41 Lax-Milgram 
theorem, Tartar’s method of oscillating test functions (see [71 pp.103] or |28j, pp.31]). We refer the reader 
to m for more details on the Tartar’s method. We also mention that all the results in Lemma 12.41 and 
Theorem 12.51 still hold for C* and B*. 

Lemma 2.7 (Cacciopolli’s inequality). Suppose that A satisfies (11.11) . and other coefficients satisfy (1 1.3 j) . 
Assume that u e £ -H/ oc (12; M m ) is a weak solution to C £ (u £ ) = div(f) + F in 12, where f = (/“) £ L 2 (12;M md ) 
and F £ L 9 (12;M m ) with q = Then for any B C 2 B C 12, we have the uniform estimate 



1 

2 




\f\ 2 dx ) 2 + rC 


2 B 


\F\ q d x y 


(2.9) 


where C depends only on fi. n, A, m, d. 

Proof. The proof is standard, and we provide a proof for the sake of completeness. Let 4> £ Cq (II) be a 
cut-off function satisfying <fi = 1 in B, 4> = 0 outside 2 B, and |V0| < 2/r. Then let tp = 4> 2 u £ be a test 
function, it follows that 


[ [AfVu% + + 2 [ [AfiN'4 + V^u?) Vfiufifdx + [ BfVu^ufcfffix 

Jn Jn Jn 

+ [ cfu^uffi 2 + \\u £ \ 2 f 2 dx = [ F a uf(j) 2 dx — [ f a S7uf(j) 2 dx - 2 [ ffiVfiuffidx, in 12. 

Jn c Jn Jn Jn 

By using the ellipticity condition and Young’s inequality with 5, we have 


^ f 4> 2 \S7u £ \ 2 dx + (A — C') f <p 2 \u £ \ 2 dx 
4 Jn Jn 

< C f \V(j)\~\u e \ 2 dx + C f 4> 2 \f\ 2 dx + f 4> 2 \F\\u £ \dx, (2.10) 

J n J £1 J £1 

where C' = C'(p, k, m, d). This together with 

j^4> 2 \F\\u £ \dx < [J (fi\u £ \) 2 dx^j 2 (^J (cj)\F\ydx S J q < \S7 {fiu £ )\ 2 dx^j 2 (J (<j>\F\) q dx} q 

< ^ / \Vu £ \ 2 f 2 dx +£ [ \V<p\ 2 \u £ \ 2 dx + C[ I \<f>F\ q dx) 2/q 

o Jn o Jn Jn 

gives (12.91) . where 2* = 2 d/(d — 2), and we use Holder’s inequality, Sobolev’s inequality, and Young’s 
inequality in order. □ 

Remark 2.8. In fact, (|2.9I) is the interior IT 1,2 estimate. By the same argument, we can also derive the 
near boundary Cacciopolli’s inequality for the weak solution to C £ (u E ) = div(/) + F in L2(4r) and u £ = 0 
on A(4r), 




1 

2 




\f~f\ 2 d x y 



( 2 . 11 ) 


We point out that the constant C in ([2.91) or (12.111) depends only on /r, k, m, d, whenever A > Aq > C'. 


Remark 2.9. Suppose that A £ A (//, t, ft), and V satisfies (11.21) and (11.41) . In view of the interior Schauder 
estimate (see [21]), we obtain 


max {HxaIIl-yy l|Vxfc||L°°(Y), [Vxdco^t)} < C(p,T, n,m,d), 


( 2 . 12 ) 
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where [Vxfc]c 0 - T (y) is the r th -H61der seminorm of (see [151 pp.254] for the definition). If A £ VMO(M rf ), 
and V satisfies ( 11 . 21 ) and (11.31) . then for any 2 < p < oo and ? £ (0,1), 

™??A\\ x k\\L°°(Y), [Xfc]c°.«(y ) 5 l|Vxfc||iP(y)} < C(n,u,K,p,<;,m,d) (2.13) 

u \ k ^ a 

can be derived from the interior W 1,p estimate (see PHHZlEI]). Note that in the case of m = 1, due to the 
De Giorgi-Nash-Moser theorem, the conditions m and m are sufficient to derive Holder estimates for 
some £ (0,1), but insufficient to get W 1,p estimates (see [T0ll221l25j ). In another special case of d = 2, K.O. 
Widman [45] obtained the Holder estimate by the hole filling technique without any regularity assumption 
on A (or see mi)- 

To handle the convergence rates, some auxiliary functions and their estimates are necessary. Let 
1 < < d, 1 < ct, /?, 7 < m, and y = x/e. Define 

i-S’fe) = »“ 7 - </(!/) - }. uriv) = vr - VT(y) - „“!>(,) (2.14) 

Note that C.E. Kenig, F. Lin, and Z. Shen [291130] showed that there exists E^A g Hp er (Y), such that 
b°f = E Hj = ~ E ilj and \\ E Uj\\L°°(Y) < C(n, u, m, d). Here we obtain similar results for Lf 7 : 

there exists F^A £ H ] jipr (Y) such that 

vr = ^ {F £ 7} ’ “ d l|rS 7 lll»(K)<Cfe W ,«,m,<i). (2.15) 

We give a proof of (12.151) for the sake of completeness. In view of (12.11) and (12.31) . we have f y U°' y (y)dy = 0 
and ^-{H“ 7 } = 0. Then there exists a unique solution Of' 1 £ H} }er (Y) satisfying 

A 0“ 7 = t /“ 7 in M d , J 0“ 7 = 0 

(see mm)- Let Ki = “ ^r(C 7 }’ and obviously F “ 7 = We mention that 6 »f 7 £ H 2 er (Y), 

which implies F^A g Hp er (Y). Next, we verify F^A is bounded, which is equivalent to 

ma?.{ll v C 7 llL°°(y)} < C{n,u,K,m,d). (2.16) 


Observe that (|2.13p gives f7“ 7 £ L p er (Y), then the estimate (12.161) follows from the L p estimates and the 
Sobolev embedding theorem, provided p > d. Moreover, we have 


d 

dyk 


{Ki} = a{C 7 ) 


d 2 

dykdyi 


W] 


ur- 


Note that ^{[/“ 7 } = 0 in M d , which implies = 0. In view of Liouvill’s theorem (see [13), we 

have -{0k 1 } = C > therefore = ^-^{C 7 } = °> and we complete the proof of (12.151) . 

In addition, we define the auxiliary functions $“ 7 and £" 7 as follows: 


At)” 7 = W? := Bg - BP(y) - B f(y)A-{x^} in R“, 

&yj 

AC° 7 = := ^ 7 - c“ 7 (y) - Bf(y)^-{xp} in M d , 


^ d^dy = 0 , 
[ ( a ^dy = 0 . 


(2.17) 


It follows from (| 2.3 [) that J Y W t a ' r (y)dy = 0 and J Y Z ai (y)dy = 0, which implies the existence of f?“ 7 and £“ 7 . 
By the same argument, it follows from (12.131) that max{||Vi?“ ||L°°(y). ||VC“ 7 ||L°°(y)} < C(fj,,u, K,d,m). 














9 


We end this remark by a summary. 
ED - El, then we have 


Suppose that A £ VMO(M d ), and the coefficients of C e satisfy 


max 


\K]\\L~(Y), linTl U~(y). l|V< 7 ||L-(y), l|VC Q7 ||L-(y)} < C{p,u,K,m,d). 


(2.18) 


In the special case of m = 1 or d = 2, the estimate (12.181) still holds without any regularity assumption on 
A. 


We now introduce the Lipschitz estimate and Schauder estimate that will be frequently employed later. 
Let L(u) = —div(AlVu) and C(u) = —div(£lVu + Vu ) + BVu + (c + A I)u. Then we have the following 
results: 

Lemma 2.10. Let Q be a bounded C 1 ’ T domain. Suppose A satisfies ED and El- Let u be the weak 
solution to L{u) = div(/) + F in LI and u = 0 on dLl, where f £ C 0,cr (Ll] M mrf ) with o £ (0, r], and 
F £ L p (D;M m ) with p > d. Then we have: 

(i) the Schauder estimate 

[Vu] c0 , CT(n) < C'{||/|| c o, CT ( Q ) + UTEp^)}; (2.19) 

(ii) the Lipschitz estimate 

l|Vii|U=o (n) < C{[f]c°- a (o) + ll-^l|Lp(f2)}) (2.20) 

where C depends on fi,r, K,,a,m,d,p and 11. Moreover, ifu = g on dLl with g £ C 1,<T (dLl\ M m ), then we have 

[Vu] c0jCT ( n ) < C{||/|| c o, CT(n ) + ||-P||LP(n) + ll5 , llci. CT (9n)}- (2-21) 

Proof. The results are standard, and we provide a proof for the sake of completeness. For (i), we refer 
the reader to ED pp.75-95] for the details. For (ii), due to the properties of Green function (denoted 
by G(x,y )) associated with L: \G(x,y)\ < C\x — y\ 2 ~ d , \V x G(x,y)\ < C\x — y\ 1 ~ d and | V x 'V y G(x, y)\ < 
C\x—y\~ d (the existence of G(x, y) is included in [26] Theorem 4.1]), we have u(x ) = — f (} V y G(x, y)f{y)dy+ 
f n G(x, y)F(y)dy. Differentiating both sides with respect to x gives 

Vrt(x) = - [ \ / x \/ y G(x,y)[f(y) - f{x)\dy - f(x) [ V x S/ y G{x, y)dy + [ V x G(x,y)F(y)dy 
Jo. Jo Jo 

= - I V x VyG(x,y)[f(y) - f(x)]dy - f(x) [ n(y)V x G(x, y)dS{y) + [ \7 x G(x,y)F(y)dy 

= ~ f V x \7 y G{x,y)[f(y) - f(x)]dy + f V x G{x,y)F(y)dy 
J o Jo 

for any x £ LI. Note that we use the integration by parts in the second equality, and the fact of V x G(x, •) = 0 
on dLl in the last equality. This implies (12.20[) . Finally we can use the extension technique to obtain 
the estimate (12.211) . (Its proof is similar to that in the proof of Theorem 11.11 and we refer the reader 
to m pp.136-138] for the extension lemmas.) □ 

Remark 2.11. Set U(x,r) = LI D B(x,r) for any x £ LI. Let ip £ Cq°(2B) be a cut-off function satisfying 
(p = 1 in B, = 0 outside 3/2 B, and |V<^| < C/r. Let w = wp, where u is given in Lemma 12.101 Then we 
have 

L(w) = div(/<^) — / • \7(p + Ftp — div(7lVv?tt) — AVuVip in Ll , w = 0 on dLl. 

We apply the estimate ()2.20l) to the above equation with r = 1, and obtain 

IIVull^^D) < C'{||u|| C -0,<7( [ /( a . i 2)) + \\u\\lp(U(x,2)) + II/IIl°°({/(x, 2)) + [/]c°’ CT (i7(a:,2)) + II^IIz,p (C/(ai,2))} 

< C{\\u\\ w i, s (u(x, 2 )) + ll/llL°°(t/(x,2)) + [f]c°’°(y(x, 2 )) + \\F\\lp(U(x, 2 ))} 


( 2 . 22 ) 
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where s = max{p, [y^] + 1} ([yzyy] is the integer part of y^yy), and we employ the fact of \\u\\c°^(u(x, 2 )) — 
C\\u\[ W i, s njt X 2 ))- On account of 

IIV^II LP(U{x,2)) < C{\\S7u\\ L Z(u( X) 4,)) + \\f\\LP(U(x,A)) + \\F\\li(U(xA))}’ 

where q = (see [21, Theorem 7.2]), we have 

ll“llw 1 - s (C/(a:,2)) < C'{||«lllV 1 - 2 (!7(a:,4)) + II/II^(£7(^,4)) + II F \\Li(U(x,A))} (2.23) 

where we use ||u \\l s (u(x, 2 )) < C{\\Vu\\ L 3 (U(x, 2 )) + \\ u \\l 2 {U{ x , 2 )) } in the above inequality (see (|3.5D ). Combining 
(12.221) and (12.231) . we have 

Il^ u llz, 00 (f7(x, 1 )) < C{\\Vu\\ L2mxA)) + \\u\\ L 2(ij{x£)) + ll/lli,°°(C/(a:,4)) + [f]c°^(U(x,4,)) + \\F\\lp(U(x,4)) } • 

Note that if U ( x , 4) C fi, then v = u — u is still a solution to L(u) = div(/) + F in 11, where u = ^ r( - x ^ udy. 
In this case, the above estimate becomes 


llVully/xq^!)) < C{||V«|| l2([/(:Ci 4)) + ||« “ ^\\L 2 (U{x,4)) + \\f \\L°°(U(x,4)) + {U{x,i)) + \\F\\lp(U(x,4)) } ■ 

If \U(x,A) n<9fi| > kq > 0, then we have IMI l 2 (u(x, 4 )) < C\\^u\\l 2 (u(xA)) due to u = 0 on dQ. (It is another 
type of Poincare’s inequality, which can be derived from the trace theorem coupled with Rellich’s theorem 
by a contradiction argument, and the proof is left to the reader.) In all, we are able to apply Poincare’s 
inequality to both of the two cases, and then obtain 


II Vt/||i,oo ([/(a:,!)) < C{ ||Vu|| L 2 ( r/ ( 3 . ! 4 )) + ||/||L°°(f/(x,4)) + [f]c 0 ^ (U(x,4)) + \\ F \\Lp (U (x ,4))}' (2.24) 

Next, we let v(y) = u(ry), where y £ U(x, 4). Then we have L(v) = div(/) + F in U(x, 4) and v = 0 on 
fl U(x, 4), where L = (ry)-^j-, f(y) = rf(ry ) and F(y) = r 2 F{ry). It follows from (|2.24l) that 

\\^ v \\l™(U(x, 1 )) < C'{||Vu||L 2 (cr(x,4)) + II/IIl°°([/(x,4)) + [f]c°^(U(xA)) + \\F\\lp(U(x,4))} ’ 


and by change of variable, we have 
l|V«|| L » W x I r )) <C , {(£ 4 |Vu| 2 dy) 


1/2 


+ 


' U (#,4 r) 

By a covering technique (shown in the proof of Theorem I4.4|) , we have 

T/2 


L°°(U(xAr)) + r<7 [f]c°’-(U(x, 4r)) +r (£ (t \ F \ P dy ) 


1/p 


\\Vu\\ L oo {u) <cUj^\Vu\ 2 dy) 


+ 


L°°(2U) 


+ r [f]c°’ CT (2U) 


2 U 


1/p 


(2.25) 


Here U is the abbreviation of U(x,r) and 2 U = U(x,2r). We mention that all the above proof is so- 
called localization argument, which gives a way to obtain “local estimates” (such as interior estimates and 
boundary estimates) from corresponding “global estimates”. The main point is based on cut-off function 
coupled with rescaling technique. So, on account of the estimate (12.19|) . following the same procedure as 
before, it is not hard to derive 


[V«W (I0 < Cr ~ a {{f 2U l V ^l 2 dy) 1/2 + \\f\\L°o {2U ) +r a [f] c0 „ {2U) +r(£ v |F| p dy) 1/P }, (2.26) 

and by (12.211) . 

[Vu] c °’°(D(r)) — \Vu\~dyj +r 1 ||fl , llL°°(A(2r)) + l|Vfif||z / oo( A( -2 r .))-|-r <T [V//] c0iCT ^ A( , 2r ,^| (2.27) 

holds for u satisfying L{u ) = 0 in D(2r ) and u = g on A(2r). We note that the estimate (12.251) is of help to 
arrive at (I2.26|) . and the extension technique (see [22] pp.136]) is used in (12.271) . The details of the proof are 
omitted. Finally we remark that ([2.271) is exactly the Schauder estimate at boundary, which can be directly 
proved (see EH Theorem 5.21]). 



















11 


Lemma 2.12. Let LI be a bounded C 1,T domain, and a € (0, r]. Suppose that A,V satisfy (11.11) and (|1.4j) . 
and B, c satisfy «• Let u be the weak solution to C{u) = div(/) + F in Ll and u = 0 on dLl, where f,F 
satisfy the same conditions as in Lemma [2.10L Then we have 
(i) the Lipschitz estimate 


l|Vu|| L ~ ( n) < C{ll/|lc°- CT (n) + 11^11^(0)}; (2.28) 

(ii) the Schauder estimate 

[Vu] c o,a^ < C{ll/llc°^(n) + ||^llLP ( n)}, (2.29) 

where C depends on /i,t, n,m,d,a,p and LI. 

Proof. The results are classical, and we offer a sketch of the proof. First we rewrite C(u) = div(/) + F as 
L(u) = div(/ + Vu) — BVu — (c + A I)u + F in Ll. It follows from (|2.20|) that 

l|Vu|| L oo (f2) < C{[f] c o, rq) + ll-F||z>(n) + IMIc°-Rn) + IMIwMfi)} 

< C{[f]c°’<r(n) + ll^lli,p(n) + 2||Vu||^ (Q )|M|^) + (2.30) 

< C{[f]c 0 >°(fi) + II-^IIlr^) + IMIl°°(£T) + IMIwurq)} + ^l|V u llL°°(f2)> 

where we use Young’s inequality in the last inequality. By the Sobolev embedding theorem we have 
IMIl°°(Q) < C|l w llw 1 >p(f 2 ) f° r P > d. This together with (I2.30P leads to 

I|Vu||l°°(q) < C{[f] c o,<r(ty + ||F|| LP(n) + ||u|| M /i. P (Q)} 

< C{\\f\\c°^(n) + ll-^IUp(n)}) 

where we use the W 1,p estimate with 1 < p < oo in the last inequality, which can be derived by a similar 
argument as in the proof of Theorem ll.il or see 11.21.221. 

It remains to show (ii). In view of (12.191) and (12.281) . we obtain 

[Vu] c0jCT( . o ^ < C'{||/||c°. CT (n) + ||F|| LP (q) + |M| L oo(n) + ||Vu|| L< x>(n) + IMlwMn)} 

< C'{ll/llc°.Rfi) + ll-FlUp(fi)}, 

where we also use W 1,p estimate in the last inequality. The proof is completed. □ 

Remark 2.13. Let u be given in Lemma 12.121 Applying the localization argument (see Remark 12.111) to 
the estimates (12.28[) and (|2.29p . we can similarly give the corresponding local estimates: 

IIVm||loo(j/) < C-jV \u\ 2 dy'j + ||/||l°°(2!7) + r<T [f]c°’°( 2 U) r {J~ 1^1 P dy S ) (2-31) 

and 

[Vu] co ,„ { u)< Cr -°[r-\fju\ 2 dy) 1 + \\f\\ L °°<2U) + r* [f] c °.° (2U) + r (j- \ F \ Pd v) ( 2 -32) 

where C depends on p,r, n,m,d,p,cr and Mq. We mention that in the proof of (I2.3ip . we also need W 1,p 
estimate like ED Theorem 7.2] for C. It can be established by using the bootstrap method which is exactly 
shown in the proof of Theorem 13.31 so we do not repeat them. The remainder of the argument is analogous 
to that in Remark 12.Ill 
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3 W 1,p estimates &; Holder estimates 

Lemma 3.1. Let 2 < p < oo. Suppose that A e VMO(M d ) satisfies (11.11) and m- Assume f = (/“) 6 
L p (12;M mrf ), and F G L 9 (12;M m ) with q = T/ien 2/ie weak solution to L £ {u £ ) = div(/) + l 7 in 12 and 

u £ = 0 on <912 satisfies the uniform estimate 

l|V« e ||^ ( n) < C'{||/|| L p ( n) + \\F\\ Lq{a) }, (3.1) 

where C depends only on p,u, n,p, q, d,m and 11 . 

Remark 3.2. The estimate (13.11) actually holds for 1 < p < oo, where q = if p > and q > 1 if 
1 < p < grr- In the case of F = o, m is shown in [6j. If F ^ 0, we can derive the above result by the 
duality argument applied in Lemma 13.71 The same method may be found in [161131] . Besides, we refer the 
reader to [317 for the sharp range of p’s on Lipschitz domains. 

Theorem 3.3. (Interior W 1,p estimates ) . Let 2 < p < oo. Suppose that A G VMO(M d ) satisfies (II.ip , (11.21) , 
and other coefficients satisfy a. Assume that u £ G H/ oc (12; M m ) is a weak solution to C £ (u £ ) = div(/) + F 
in 12, where f G L p (12;M mrf ) and F G L' ? (12;M m ) with q = Then, we have |Vit e | G L P (B) and the 

uniform estimate 



for any B C 2 B C 12 with 0 < r < 1, where C depends only on p , u, k, A,p, m, d. 

Proof. By rescaling we may assume r = 1. In the case of p = 2, (13.21) follows from Lemma 12.71 Next, we 
will prove (13.21) for p G [2,p*], where p* = d ^2k 0 ’ anc ^ — &0 + 1- To do so, let w £ = pu £ , where 

p G C'o°(2R) is a cut-off function satisfying tp = 1 in B, p = 0 outside 3/2 B, and |V<p| < C. We rewrite the 
original systems as 

—div(yl £ Vi(; £ ) = div(/<p) — / • V<p + Fp + F in 12, 

where 

F a = divfl rfwl - Afv<pu%) - Af VpS/4 - VfVpuV - Bfv^p - cf - A 
Hence it follows from (12.91) . (13.11) . (13.61) . and Holder’s inequality that 

||Vu £ || lp(b) < C{\\u £ \\ LP{2B) + \\Vu £ \\ L1(2B) + ||^e||L9(2B) + ||/||lp(2B) + \\F\\li(2B) } 

< C {\\y U £ \\ L q(2B) + ll^e I!l 2 (2B) + II/IIlp(2S) + I|T||l9(2B) } 

< C | ||^e|lL 2 (2 , =o+ 1 B) + 11/ II Lr{2 k o+ 1 B) + \\F\\L,i{2 k o+ 1 B) | ■ 

We first check the special case of p = p* to obtain the final step k$ of iteration, and then verify the above 
inequality for any p G [2, Second, it is not hard to extend the range of p’s to [2,p*] by at most ko times 
of iteration. The rest of the proof is to extend the p’s range to p* < p < oo. Indeed it is true, since 

l|V« £ ||i>(£) < C'{||Vu e || W(2fl) + ||^e|| l 2 (2B) + II/IIlp(2B) + I|T’||l9(2B) } 

and q < d G [2,p*], which is exactly the start point for iterations due to the previous case. 

Hence, let N = A k ° , we have proved 

\\^ u e\\lp(B) < C{ll' u e||L 2 (iVB) + ||/||lp(JVB) + Ill’ll Ls(NB)} (3-3) 

for any 2 < p < oo in the case of r = 1. We remark that (i) the estimate (13.31) uniformly holds for e > 0; (ii) 
the constant in (|3.3I) can be given by C < C(p, u, m, d,p){ ||^4.||z. oo (M rf ) + ll^llL°°(R d ) + ll-®llL°°(R d ) + ll c llL oo (R d ) + 
A} fc o +2 . The two points make the rescaling argument valid when we study the estimate (13.21) for 0 < r < 1. 
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We now let v £ (x) = u e (rx), where x £ B\ and r £ (0,1). Hence we have 
Cl(v £ ) = — div[H(rx/e)Vt! £ + V(rx/s)v £ ] + B(rx/e)Vv £ + c(rx/e)v £ + Xv £ = div(/) + F in NB\, (3.4) 
where 


V(x) = rV(x), B = rB(x), c = r 2 c(x), X = r 2 X, f = rf(rx), F = r 2 F(rx). 

It is clear to see that the coefficients of C satisfy the same assumptions as C in this theorem. Set e' = e/r, 
and applying (13.311 directly, we obtain 

IIVMMffi) — ^ {\\ v re'\\L 2 {NB 1 ) + H/llLP(iVBi) + 1111 (JVBi) } > 
where C is the same constant as in (13.311 . This implies 

||Vu e || iP ( Br ) < C{r 1+ p 2 \\u e \\L 2 (NB r ) + l|/IUp(JVB r ) + r 1+ p q \\F\\ L q( NBr )}- 
Finally, for any B with 0 < r < 1, we choose the small ball with r/N radius to cover B r . Hence we have 

\\^ u e\\LP{B r ) < C{r 1+ p 2 ||M e || L 2( B2r ) + ||/|| LP ( B2r .) + r 1+ P q \\F\\Lt(B 2r )}i 

and this gives (13.211 . We complete the proof. □ 

Remark 3.4. Here we introduce two elementary interpolation inequalities used in the above proof. Let 
u £ W 1,p (Ll] M m ) with 2 < p < oo, then for any 5 > 0, there exists a constant Cs depending on 5,p, d , m and 
Ll, such that 

IMIlp(Q) < ^l|Vn|| iP (Q) + Cs\\u\\ L 2^. (3.5) 

The estimate (13.511 can be easily derived by contradiction argument (or see [lj). As a result, we have 

IMIlp(O) < CIMliW'h^) — C , |I^ 7 u IIl5(o) + f 1 I^||l2(q) (3.6) 

for 1 < p < oo and q = where C depends on p, d , m and Ll. 

Corollary 3.5. Suppose that the coefficients of C £ satisfy the same conditions as in Theorem 13.31 Let p > d 
and cr = l — d/p. Assume that u £ £ H/ ) (Ll-,W n ) is a weak solution of C £ {u £ ) = div(/) + F in Ll, where 
f £ L p (Li]M md ) and F £ L q (Ll]R m ) with q = ^ > \. Then we have 

W u e\\c°’°(B) < Cr~ a { (J-^ \u £ \ 2 dx^j 2 +r(£^\f\ p dxY + r 2 ( K j-^\F\ q dx s j q ^ (3.7) 

for any B C 2 B C LI with 0 < r < 1. In particular, for any s > 0, 

IMIl°°(b) < C \{i \ue\ s dxy +r(i \f\ p dx'j p + r 2 (i- \F\ q dx^ q j , (3.8) 

v ’J 2B J 213 J 213 j 

where C depends only on p,co, k, X,p,m, d. 

Proof. Assume r = 1. It follows from the Sobolev embedding theorem and Remark 13.41 that 
ll' u ellc 0 " CT (B) < C\\u £ \\ w i, P ^ < C'j|Vn e || BP ( B ) + C\\u £ \\l2(b)- 
Then it follows from (13.21) and rescaling arguments that 

\\ u e\\c q ’°{B) < Cr~ u \u e \ 2 dx\ +Crl a \(j~ 2 \f\ Pdx ) V + r (/f B ’ 
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where a = 1 — d/p. Moreover, for any x £ B we have 


u e (x) | < | u e (x) - u £ 


\Ur 


< [u £ ] 


C°’ a (B) 


r a + 


£ 


M\ 


\ 1/2 


from Holder’s inequality. This gives 

IMIl°°(B) < c (S ^ \u E \ 2 dx^ + Cr | (j- B \f\ P dx^j P + r (J~ b ? | • 

Moreover, by the iteration method (see [2Q pp.184]), we have (13.81) . □ 

Lemma 3.6. Suppose that the coefficients of C £ satisfy the same assumptions as in Theorem 13.31 Let 
1 < p < oo, f = (/“) £ L p (Q-,M. md ). Then there exists a unique u £ £ B / 0 1 ’ p (S7;R m ) such that C e {u £ ) = div(/) 
in and u £ = 0 on dLl, whenever A > Aq. Moreover, the solution satisfies the uniform estimate 


l|Vu E ||LP(n) < C\\f\\ L p(Q), (3.9) 

where Ao is given in Lemma |2~H and C depends on p,co,K,X,p,d,m and Ll. 

Proof. In the case of p = 2, it follows from Theorem 12.51 that there exists a unique solution u £ £ 17g(ll;M m ) 
satisfying the uniform estimate || ||z/ 2 (r 2 ) < C\\f\\L 2 (n)- For P > 2, the uniqueness and existence of the 

weak solution is reduced to the case of p = 2. We rewrite the original systems as 


L £ (u £ ) = div(/ + V £ u £ ) - B £ Vu £ - (c £ + A)it e . 

Applying (12.81) . (13.11) and Sobolev’s inequality, we obtain 

||Vn £ || L Pfe 0 (n) ^ C{II/IIl p,c o( si) + ll' u ellz, Pfc o(n) + ll^ n £llL Pfe o- 1 (n)} — C{\\f\\ L Pk 0 ^ + ll^ 7 ' u ellz p fco- 1 (n)} 

< C , {||/llL Pfc o(f 2 ) + ll/lli 2 (n)} < C|l/llL Pfc o(f 2 )’ (3.10) 

where pk 0 = d-2ko 1 anc ^ * s a P os ibi ve integer such that fco < d/2 < ko + 1. We claim that (I3.10P holds 
for any p € [2,p^ 0 ]. Indeed, let T(f) = Vu £ , then together with ||T , || i 2_ ) . i 2 < C and \\T\\ L Pk 0 _^. L Pk 0 < C, 

the Marcinkiewicz interpolation theorem gives ||T||< C (see SU)- Moreover, for any p > pk 0 , we still 

have ||Vn £ || L p (f2) < C{\\f\\ LP ( n) + || Vu £ || L9(n) } < C||/|| LP(n) , since q < d. 

By the duality argument, we can derive (13.91) for p £ (1,2). Let h = (/if) £ Cg (LI; R md ), and v £ be 

the weak solution to C*{v e ) = div(/i) in 12 and v £ = 0 on <912. Hence, in view of the previous result, we 

have ||Vu £ || L j/(q) < C\\h\\ LP i^ for any p' > 2. Moreover, if / £ Cg (12; M md ), there exists the weak solution 
u £ £ i7g(12;M m ) to the original systems. Then it follows from Remark 12.31 that 

/ Vu £ hdx = — u £ C*(v £ )dx = — C £ (u £ )v £ dx = / fWv £ dx. 

Jn Jn Jn Jn 

This gives ||Vit £ ||x /P (Q) < C'H/Hlp^), where p = p'/(p' — 1). By the density argument, we can verify the 
existence of solutions in Wq’ p (LI) for general / £ L p (Ll;W nd ), as well as the uniqueness for 1 < p < 2. The 
proof is complete. □ 

Lemma 3.7. Suppose that the coefficients of C s satisfy the same conditions as in Lemma [3.61 Let 1 < p < 
oo. Then for any F £ L q (Ll; M m ), where q = pdf (p + d) if p > d/(d— 1), and q > 1 if 1 < p < d/(d—T), there 
exists a unique solution u £ £ Wg’ p (12; M m ) to C £ (u £ ) = F in LI and u £ = 0 on <912, satisfying the uniform 
estimate 

l|Vu £ ||iP(f2) < C\\F\\ L q(n), 

where C depends only on p,uj, n, X,p,q,d,m and 12. 


(3.11) 
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Proof. We prove this lemma by the duality argument. The uniqueness is clearly contained in Lemma 13.61 
and the existence of the solution u £ follows from the density and Theorem 12.51 The rest of the proof is to 
establish (13.111) . 

Consider the dual problem for any / € Cq^M" 1 ^), there exists the unique v £ £ iLo(12;M m ) to £*(u £ ) = 
div(/) in 12 and v £ = 0 on d£l. Note that it follows from Lemma iThl that ||Vu e || iP '(n) — C\\f\\ L p'(ay Then 
we have 


Vu £ fdx = - u £ C*(v £ )dx = - C £ (u £ )v £ dx = - Fv e dx , 
! in Jn in 


and 


\7u e fdx 


— Il^llw(n)||u c ||^ (n ) — ||Vw £ || LP '^ < CH-FH^^II/ll^y^). 


Note that ^ = jy — ^ if p' < d, 1 < q' < oo if p' = d, and q' = oo if p' > d. In other words, q = -p^ if 
P > 3 ^ 1 , and < 7 >lifl<p< Finally we obtain ||Vu £ || L p(o) < C'||F|| L9 (n). □ 


Proof of Theorem 11.11 In the case of g = 0, we write v £ = u £} i + u £: 2 , where u £ q and u £ 2 are the 
solutions in Lemma 13.61 and 13.71 respectively. Then we have 


I|Vu £ || LP (q) < ||Vu £i i|| iP (Q) + ||Vu £i2 ||i,p(q) < C'{||/||z,p(q) + ||i ? || L9 ( r2 )}. (3.12) 


For g / 0, consider the homogeneous Dirichlet problem C £ {w £ ) = (1 in Q and w £ = g on dLl, where 
g € B 1 ~ 1 / p ’ p (dLl; M m ). By the properties of boundary Besov space, there exists G £ W 1,p (H;M m ) such that 
G = g on dQ and ||G , ||u/ 1 '? , (n) — ^'\\9\\b 1 ~ 1 /P’P(dii)- Let h £ = w £ — G, we have 

C £ {h £ ) = div(7L £ VG + V e G) - B £ VG - (c £ + A )G in fi, h £ = 0 on dQ. 

Recall the case of g = 0, in which there exists the unique weak solution h £ £ VF 0 1,p (n; M m ), satisfying the 
uniform estimate ||Vh £ ||/>(o) < C'||G|| l yi,p(Q) < C\\g\\ B i-i/ p , p / d ^ for 1 < p < 00 . This implies 

l|Vu; £ || z , P (f 2 ) < ||Vh £ || i p(Q) + ||VG|| L p(n) < C\\g\\ B i-i/ p , p ^ an y (3.13) 

Finally, let u £ = v £ + w £ . Combining (13.121) and (|3.13l) . we have 

l|Vu £ ||z,p(n) < C{\\f\\ LP (ty + ||T’||z,9(f2) + IIpIIbi-i/p.p(9Q)}) 

where C depends only on //, w, k, A ,p, q, d, m and Q. We complete the proof. □ 

Corollary 3.8. Suppose that the coefficients of C £ satisfy the same conditions as in Theorem II .11 Set d < 
p < 00 and a = 1 - d/p. Let f = (/“) £ L p (D;M md ), F £ L9(D;R m ) with q = and g £ C 0 ’ 1 ^; M m ). 
Then the unique solution u £ to (EZD satisfies the uniform estimate 

ll w e||c 0 ’Rf2) < C{\\f\\ LP(n) + ||T||l 9(Q) + Ibllc'O.Ran)}, (3-14) 

where C depends only on p,,uj,n,\,p,q,d,m, and f i. 

Proof. Due to the extension theorem (see ] 22| pp. 136]), there exists an extension function G £ C' 0,1 (D;M m ) 
such that G = g on dfl and ||G|| c o.i ( q) < C\\g\\ c o,i^y This also implies ||G|| w i.p( n) < C\\g\\ c o,i( m) for any 
p > 1. Let v £ ,w £ be the weak solutions to the following Dirichlet problems: 


(i) 


C £ (v £ ) = div(/) + F in 0, 

v £ = 0 on dLl, 


(ii) 


C £ {w £ ) = 0 

w £ = g 


in 12, 
on dLl, 


(3.15) 


(fi) - 


respectively. For (i), it follows from the Sobolev embedding theorem and Theorem 11.11 that ||u £ ||co,. 
C|I v ^IIlp(o) < C{\\f\\ L p(n) + ||-F||z,9(n)}- For (ii), by setting h £ = w E - G, we have C £ {h £ ) = -C e {G) in 
Q and h £ = 0 on dQ. Therefore ||/i £ ||c-o.*(n) < C'||G|| w i,p( n ) < C\\g\\ c o,i {an) , which implies ||w £ ||c°.<qn) < 
C , ||5 f llc°. 1 (9n)- Let u £ = v £ + w £ . Combining the estimates related to v £ and w £ , we derive the estimate 

(l3~1~4l) . □ 














16 


Remark 3.9. Assume the same conditions as in Corollary 13.81 let u £ be a weak solution to (II.7|> . Then by 
the same argument as in the proof of Theorem 13.31 we obtain the near boundary Holder estimate 

\\Ue\\co^(D(r)) < Cr~ a | (^ {2 ^ \u E \ 2 d X y + r llidlco.i (A(2r)) + r ^ \f\ P dxj P + ^ ) \ F \ qdx ) * | > 

(3.16) 

where cr = 1 — d/p, and C depends on y, ui, k, A ,p, d, m and fl. 

Remark 3.10. In the following, we frequently use the abbreviated writing like f^F(-,y) = f^F(x,y)dx. 

Theorem 3.11. Suppose the same conditions on C £ as in Theorem, II.1L and A > Ao- Then there exists a 
unique Green’s matrix G £ : 12 x Q —>• M ™ 2 U {oo}, such that Q e (-,y ) G H 1 (fl \ B r (y)-,W n2 ) n W 0 1 ,s (f2; M™ 2 ) 
with s G [1, jry) for each y G fl and r > 0, and 

B :[Gf(;y)A\ =^(y), 

Particularly, for any F G L 9 (J7;M m ) with q > d/2, 

u/(y)= [ Qf a {x, y)F a (x)dx 

Jci 


(3.17) 


(3.18) 


satisfies C £ (u £ ) = F in Gt and u £ = 0 on <912. Moreover, let *Q e {-,x) be the adjoint Green’s matrix of G £ (-,y), 
then G £ (x, y) = [ *Ge(y, x) ]* and for any a, a' G (0,1), the following estimates 


\Gs(x,y)\ < 


C 


\ x - y I 


d—2 


mm 


in jl, 


d° 


d°, 


dfdf 


\x — y\ a ’ \x — y\ a ' ’ |a; — y\ a+ ' 




(3.19) 


hold for any x,y G Gl and x / y, where d x = dist(x, <912), and C depends only on y, oj, n, A, d, m and 12. 

Lemma 3.12 (Approximating Green’s matrix). Assume the same conditions as in Theorem 13.111 Define 
the approximating Green’s matrix G p>e (-,y) as 


B l\Gl, e {;v)A = / u^dx, Vug H'(Gl-,R m ), 


o(y) 

where 1 < 7 < m, and Gl p (y) = 12 n B p (y). Then if \x — y\ < d y /2, we have the uniform estimate 

C 


I Gl £ {x,y)\< 


\x - y\ 


d—2 ' 


v P < \x-y\/4, 


where C depends only on y,tw, k, A, d,m and 12. Moreover, for any s G [1, -jxj), we have 

SU P ll G 'p,e('’y)|lvy 1 ’ s (m - C(y,u,K,\,d,m,s,n,d y ). 

p >0 0 ' 


(3.20) 


(3.21) 


(3.22) 


Proof. First of all, we show G Pt£ (x,y ) = [G/,® £ (x,y)] is well defined. Let I(u) = ^ u 1 dx , then I G 

F[~ 1 (Gl;M m ) and |I(u)| < C\Gl p (y)\^ 1 / 2 * \\u\\ H i^ with 2* = ^4j. It follows from Theorem 12.51 that there 
exists a unique Gfp %e (-,y) G i7d(fi;M m ) satisfying (|3.20D and 

l|VGp ie (-,y)||z, 2 ( n ) < C\Tl p (y)\~^. (3.23) 

For any F G C//°(Gl] M m ), consider C £ (u £ ) = F in 12 and u £ = 0 on <912. There exists the unique solution 
u £ G TTq (0; M." 1 ) such that 


f FGl £ (-,y) = B £ (u £ ,Gl £ (;y)] = B*[G^(-, y),u £ ] = l 
J Cl J Cl 


U’ 


ftp(y) 


(3.24) 
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Suppose supp(F) C B C Q, where B = Br{d). Then it follows from (12.81) . (13.8j) and (13.241) that 


[ F G/ £ (-,y) < ||u e || i oo( 1 / 4B ) < C 
Jn 


for any p < R /4 and p > d/2. This implies 


u f 


i i [■ i 

| 2\2 , r >2 / 


1/2B 


+ R“ 


1/2B 


l*T 


< Ci?" 


|F| P 




<Ci? 2 - d , Vi?<d„ V 9 G [1,—). 

Now we turn to (13.211) . Set r = \x — y |, and r < d y /2. In view of (13.201) . G pe (x,y ) actually satisfies 
C e [CT p ^.,y)] =0 in Bz(y) \ Bn(y). By using (13.81) again, we obtain 


\GJ, e (x,y)\ < C-f \Gl e (;y)\<C-f |G^(-,y)| < C\x - y\ 

JBr(x) J B2r(y) 


2 -d 


for any p < \x — y|/4, where C depends only on p, u, k, A, d, m and Jl. 
Then we will prove (13.221) . Step one, we verify the following estimates, 


[ \VGl E (-,y)\ 2 <CR 2 - d , f \G'J IE (-,y)\ 2 * <CR~ d , V p > 0, V R < d y /4. (3.25) 

Jfl\B{y,R) Jn\B(y,R) 


ln\B(y,R) 


On the one hand, let G Cq (O) be a cut-off function satisfying p = 0 on B{y, R), cp = 1 outside B(y, 2R), 
and |V<p| < C/R. Choose u = cp 2 Gp t e(-,y) in (|3.20l) and A > Aq. It follows from p.101) and (13.211) that 


[ <p 2 \VG/ j£ (;y)l 2 <C f |Vp| 2 |G^(-,y)| 
Jn Jfi 


9 C 
2 < 


R 2 


'2 B\B 


- y\ 2{2 ~ d) < CR 


2—d 


V p< R/4. (3.26) 


On the other hand, it follows from (|3.23|) that 


|VG^(-,y)| 2 < / \VGl £ (-,y)\ 2 <CR 2 ~ d , Vp>R/4 


ln\B(y,R ) 


Thus we have the first inequality of (13.251) for all p > 0. 
For the second estimate in (13.251) . we observe 


WGl e {-,y)\ 2 ~ <C( I \V{pGl £ (;y))\ 2 ) d ~ 2 <C( I \ VcpG/^, y)\ 2 + |<pVG^(-, y)| 


d 

|2\ d—2 


< Ci? 


-d 


(3.27) 

for any p < R/ 4, where we use Sobolev’s inequality in the first inequality and (|3.26l) in the last inequality. 
We remark that the constant C does not involve R. In the case of p > R/A , since Gp i£ (-,y) = 0 on 50, we 
have 


/ 

do' 


|G£ £ (-,y)f < / |GX £ (-,y)| 2 * <C( / |VG2, e (-,y)| 




2 \ d-2 




< Ci? 


—d 


ln\B(y,R) 

where we use Sobolev’s inequality in the second inequality and (|3.23l) in the last inequality. This together 
with (|3.27j) leads to 


\Gl £ (-,y)f <CR~\ Vp> 0, \/R<d y / 4. 


ln\B(y,R) 


We now address ourselves to the uniform estimates of Gp, e (-, y) and VGp, e (-, y) with respect to parameter 
p. In the case of t > (d y / 4) 1_d , we obtain 


| {x G fl : |VGp e (-, y)| > f} | < CR d + f -2 / |VG£ e (-, y)| 2 < V p > 0. 

Jn\B(y,R) 


(3.28) 
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For t > ( dy/4:) 2 , it follows that 

\{xen:\Gl E (;y)\>t}\<Ct~^, V p > 0. 

Then in view of (|3.28p and (13.291) . we have 

r roo d 

/ | (P 0 J-,y )\ 8 < Cdt^-^ + C / t 3 - 1 -t-^dt < C\df^ +d s v ^ +d } 

JS 2 J(d,JiV-d 


for s G [1, 3 ^ 2 ), and 


|VG^(- ! 2 /)| s <C[df- d )+df- d ) +d ] 


(3.29) 


for s G [1, where C depends only on fj,,u, k, X,d,m,s and 12. Combining the two inequalities above, 
we have (13.221) . and the proof is complete. □ 


Proof of Theorem 13.111 From the uniform estimate (1 3.22 [) . it follows that there exist a subsequence of 
{ Gl p n ,e{^y)}n= 1 and such that for any s G (1, gry), 

G 7 ni£ (-,y) 5 e 7 (-,l/) weakly in W 0 1 ,s (12;M m ) as n-s- 00 . (3.30) 


Hence, we have 

BeW(-,y), 0 ] = linr B*[G^, e (-,y)^] = km / ^ = ^(y) 

for any (j) G Wq ,p (12; M m ) with p > d, where we use the definition of the approximating Green’s matrix. Due 
to Theorem ll.il there exists the weak solution u £ G Wq ’ p (12: M m ) satisfying C £ (u E ) = F in 12 and u £ = 0 on 

P 

<912 for any F £ L 2 (12; M m ) with p > d. Thus we obtain 

u2(y) = B* £ [g](.,y),u E ] = B £ [u £ ,g](;y)}= [ GK;y)F. 

Jq 


We now verify the uniqueness. If Gl (-, y) is another Green’s matrix of we then have uj = Jp Gl (-, y)F. 
It follows from the uniqueness of the weak solution that f n \Gl (■, y) — Gl (*, y)] F = 0 for any F G L 2 ( 12 ; M m ), 
hence Gl{-,y) = Gl{-,y) a.e. in 12 . 

Next, let *G B)£ (-,x) denote the approximating adjoint of G Pt£ (-,y), which satisfies 

B £ [*Gl E (;x),u] = -f u e , VuG lFg(12;M m ). (3.31) 

By the same argument, we can derive the existence and uniqueness of *G £ (-,x), as well as the estimates 
similar to (13.211) and (I3.22|) . Thus for any p, g > 0, we obtain 



B * £ [Gl £ (;y),*GU;x)\ 




B £ [*Gl £ (-,x),Gl e (F, 



Note that C £ [*G e g £ (-,x)\ = 0 in 12 \ B g {x) and <C*[Gp,e(‘> y)] = 0 in 12 \ B p {y). In view of Corollary 13.51 
*G d ge (-,x) and Gp >e (-,y) are locally Holder continuous. Therefore, we have *Gl e (y,x) = Gl^(x,y) as p —> 0 
and g —> 0, which implies G £ (x,y) = [*G e (y, x)}* for any x,y G 12 and x / y. 

Let r = \x — y\ and F G C'^°(12r(x)). Assume u £ is the solution of C £ {u £ ) = F in 12 and u £ = 0 on <912. 
Then we have u e (y) = J n Ge(z,y)F(z)dz. Since C £ (u £ ) = 0 in 12 \ 12r^, it follows from Corollary 13.51 that 

\u e (y)\ < c(-f \u £ {z)\ 2 dz) 2 <Cr x ~U [ \u £ (z)\ 2 *dz) 2 * 

Vfi r ( s ) ' Wo ' 
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< Cr 1 


\\7u £ (z)\ 2 dz) < Cr l s\\F 


2d 

TT(f2) 


< Cr 1 


ll^llL^n, (x)) > 


where we use Holder’s inequality in the second inequality and Sobolev’s inequality in the third inequality, 
as well as the estimate (13.111) with p = 2 in the fourth inequality. This implies 


\Ge(z,y)\ 2 dz ) 2 < Cr 


ft r (x) 


„2 -d 


(3.32) 


Note that C e [G e (-,y)] = 0 in fi \ B(y,r) for any r > 0. So in the case of r < 3 d x , it follows from (13.81) and 
(13.321) that 


\Ge(x, y)\<C[4- \G e (z,y)\ 2 dz) 2 < 


Q r (x) 


C 


\x — y\ d 2 

For r > 3 d x , in view of (13.161) and (13.321) . for any a £ (0,1), we have 


x — x\ 


\Qe{x,y)\ = | Ge(x,y) - Ge{x,y) | < C{- -—) 


\Ge{z,y)\ 2 dz) 2 < C 


Q r (x) 


\x -y\ d - 2 +°' 


where x € dGl such that d x = \x — x\. By the same argument, we can obtain similar results for *Q e (-,x). 
Since \G s (x, y)\ = \*G e (y, *)| for any x,y G fl and x / y, the following results are easily derived by the same 
arguments: 


\Gs(x,y)\ < 


CdZ 


— min 11, d ° x } 
7 l x — v a ■> 


x - y\ c 


\x — y| rf - 2 +°'' 

for any a, a' £ (0,1), where C depends only on p, w, k, A, d, m, a and Q. The proof is complete. □ 

Remark 3.13. We will see in Section 4 that the estimates (13.191) actually hold for cr = a' = 1, which are 


\G s (x,y)\ < 


C 


\x - y | 


d—2 


mm 


i 1 - 


dii 


d x d x dy 


\x-y\' \x-y[ \x-y\ 2 


Let r = \x — y |, and y £ dft such that d y = \y — y\. In the case of d y < r/6, due to G e {x, ■) = 0 on <912, we 
have 


\Ge(x,y)\ = I Ge{x,y) - Ge(x,y) | 


< \\^Ge(x,-)\\L^(nr(y))\y -y\ < — (/ \Gs(x,z)\ 2 dzY < C \^J 

v n r \Jnr(y) 1 \x-y\ d 1 


(3.33) 


where we employ the estimate (14.361) in the second inequality, and (13.321) in the last one. In the case of 
d y > r/6, we can straightforward derive the above estimate from the estimate \G £ (x,y)\ < C\x — y\ 2 ~ d - 
Similarly, we can derive 

\Ge{x,y)\ = \*Ge(y,x)\ < Cd x \x - y\ l ~ d ■ 

Then we plug the above estimate back into the last inequality of (13.331) . and obtain \G £ {x,y)\ < Cd x d y \x — 
y\~ d . 


Remark 3.14. The main idea in the proofs of Theorem 13.111 and Lemma 13.121 can be found in [3ll26]. We 
comment that the indices a and a' £ (0,1) can be equal, which actually come from the Holder estimate 
with zero boundary data. Equipped with the estimate (|3.19l) . it is possible to arrive at the sharp Holder 
estimate with nonzero boundary data. 



































20 


Proof of Theorem II.2L We first assume that u e \ satisfies C £ {u £ \) = 0 in fl and u £t i = g on dGl. Let v 
be the extension function of g, satisfying Av a = 0 in Gl and v a = g a on dGl. For any i6l!, set 
B = B(x, d x ). We have the estimate 

|Vu(x)| < C (^ b |Vv| 2 dy) 2 < \v(y)-v(x)\ 2 dyY < C'd£ -1 [u]c*,* ( n ) < Cd a x ~ l \\g\\ c o,a [dn) (3.34) 

for any a € (0,1), where we use the (interior) Lipschitz estimate (12.251) in the first inequality, Cacciopolli’s 
inequality in the second inequality, and the Holder estimate: [u]co,<t(q) < C\\g\\co,<r(Qty in the last inequality. 
By normalization we may assume ||g||c°’ CT (an) = 1. Let w e = u £j i — v. then C £ {w £ ) = —C £ {y) in Gl and w £ = 0 
on dGl. It follows from (13.181) that 


w £ (y) = — / VG £ (x,y)[A £ S/v + V e v\dx - / Q e (x, y) [B e \7v + (c £ + A)u] dx, 

which implies 

\w £ (y)\ < C [ \\/G £ (x,y)\[d x } a ~ 1 dx + C f |£ e (z,y)|[4] CT_1 dz + <7 /" (\VQ e (x, y)\ + \G e (x, y)\)dx 

JQ JQ J£l 

= : I\+ h + h, 


where we use the estimate (13.341) . 

To estimate I\, set r = d y /2. It follows from (12.91) and (13.191) that 


/ \S7 x Ge{x iy )\[d x Y- l dx < CV(2-V) d 

'B{y, r) P) 


2 2 „cr—l 


B(y,2-jr)\B(y,2-j~ 1 r) 


|V a ;Ge{x,y)\ dx) r 


< C^2(2-^r) d ~ 1 

3=0 


\Ge(x,y)\ 2 dx] r a ~ 1 < Cr a , 


B(y,2-j+lr)\B(y,2-j-lr) 


Next, we address ourselves to the integral on Gi \ B{y,r). Let Q be a cube in M. d with the property that 
3 Q C fl\{y}, and l(Q), dist(Q, dGl) are comparable, where l(Q) denotes the side length of Q (see (U] pp.167]). 
Thus, for fixed z G Q there exist ci, C 2 > 0 such that c\\z — y\ < \x — y\ < C 2 \z — y\ for any x G Q, and we 
have 



\V x G e (x,y)\ [dtf-'dx < C[l{Q)Y~ 2 \Q\ 



\Ge(x,y)\ 2 dx"j 


1 

2 


< c[im a+ai ~ 2 \Q\ 


z — y\d—2+ai+<T2 


< Cr a2 



cr+cr i —2 
d—2+ai +<72 


dx, 


where we use the estimate (12.91) in the first inequality, the estimate (13.191) in the second one, and the 
Chebyshev’s inequality in the last one. (Note that o\ and cr 2 will be given later.) By decomposing f l\B(y, r ) 
as a non-overlapping union of cubes Q (see m pp. 167-170]), we then obtain 

r \(j 10-+0-1—2 

V x Gs(x,y)\[d x r- l dx<Cr"> / , / xJ iw_ 2+a - +fT dx =: In + J 12 + J 13 . 

.7s 0 uSius 2 [r + \x-y\) d 2 +^ 1 +^ 

Note that we add the additional distance r in the denominator in the second inequality and therefore the 
corresponding domain of integral becomes the union of So = B(y, r), Si = and S 2 = U^S 0 up o Hj )m , 

where 




ft\B(y 


Qj = Qn{2 j r<\x-y\<2 j+1 r}n{2 j r<d x <2 j+1 r + 2r}, 
n jjm = nn{2 j r<\x-y\<2 j+1 r}n{2- m - 1 (2 j r)<d x <2- m (2 j r)}. 
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Then a routine computation gives rise to Jn < Cr a , 


112 

I VS 


< 


< 


Cr a * ■ 


(2 b 


icr+cr i—2 


3=0 


{23 r) d ~ 2+rJ1+a2 


• (2 ir) d < Cr a ^ j {2 j Y~ (T2 < Cr a , 


3=0 


°° 00 ro— m— 1 (<yj r \]cr+ai— 1 °° 

C'r' 72 • EE- (2i r )rf-2+o-i+a 2 • (2 ir ) <C^ 5^(2 ) + 1 


j =0 ra =0 


m=0 


J2(2 r y-' J2 < Cr a , 

3=0 


provided we choose o\ , o-i G (0,1) such that <7i + a > 1 and 02 < cr. Combining the above estimates, we 
obtain /1 < C[d y ] a . In view of (13.191) . we obtain 

h < C[i y \° j n |x _ 1 y| ,_ 1 dx < C[V, 

and /3 < C(/i + 12 ). Hence, for any y € H we have 


|^e(y)| < A + I 2 + L 3 < C[d y ] a . (3.35) 

Consider three cases: (1) \x — y\ < d x /4] (2) \x — y\ < d y / 4; (3) |x — y\ > max{c4/4, d y /4}. In the first 
case, let r = d x . In view of (|3.34l) and (|3.35l) we first have 

sup |Vu(z)| < Cr 0--1 and sup \w £ {z)\ < Cr a . (3.36) 

z&B(x,r/ 2) z£B(x,r/ 2) 

It follows from (13.71) . (13.351) and (13.361) that 

\w £ {x) - w £ (y)\ < [w £ \ C o,°( B (x idx /i))\x - y\ a 

< C\x — y\ a ^r~ a (J- \w £ \ 2 dz^j 2 + r ' lr 7 ( y j~ (|Vu| p + |u| p )dz^J P + r 2 ~ a (^J- (|Vu | 9 + \v\ q )dz^j q 

< C\x — y\ a , 

where we also use the fact ||u||^cx.(q) < C (the maximum principle) in the last inequality. It is clear to see 
that we can handle the second case in the same manner. 

In the case of (3), we derive 


\w £ {x) - w £ {y )| < |w £ (x)| + \w £ {y)\ < C\x - y\ a 

from (13.351) . Thus we have proved ||«; e ||c°> CT (fi) — L1|5 , IIc°’ ct ( 9Q)- This together with ||u||co,<t(q) < C , ||5||c°. ct (Q) 
gives 

||'«e,i||c 0 ^(a) < ll^ellco^ia) + IMIc°.°-(fi) — C\\g\\c°^{au)- 

In addition, assume that u £j 2 satisfies C £ {u £y 2 ) = div/ + F in 0 and u e> 2 = 0 on <90. It follows from 
Corollary EJ2 that ||u £ ,2||c°. ct (Q) < C'{||/||LP(fi) + H-FlU^fi)}- Let u £ = u £ g + u £ , 2 , we finally obtain ([L9|) and 
complete the proof. □ 


4 Lipschitz estimates &; Nontangential maximal function estimates 

Lemma 4.1. Suppose A € A {/j,,t,k). Let p > d and v G (0, 77 ]. Assume that f = (/“) G C 0,v {Ll] R md ), 
F G L p (0; M m ) and g G C 1,l> {dQ; M m ). Then the unique solution u £ to L £ {u £ ) = div(/) + F in 0 and u £ = g 
on <90 satisfies the uniform estimate 

I|Vu £ ||loo(q) < C'{||/|| C o,^(f 2 ) + \\FWw) + llffllci.^ian)} > (4-1) 

where C depends only on g, r, k, m, d, rj,p, v and O. 
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Proof. See m Remark 16]. In fact, (14.11) is a updated version of [3] Theorem 2]. □ 

Lemma 4.2. Suppose A £ A Let T £ (x,y ) denote the fundamental solution of L e , then we have 

max{|V x r e (x, z)\, \\7 z T £ {x,z)\}<C\x-z\ l ~ d , \V X V z V £ (x, z)\ < C\x - z\~ d , (4.2) 

where C depends only on p, r, k, m, d. 

Proof. See m pp. 6 ], □ 

Lemma 4.3. Suppose A £ A(/z, r, k). Let F £ L p (fl;M m ) with p > d, f £ C 0,fr (H; M md ) with a £ (0,1), and 
u £ £ Hl oc (fl;R m ) be the weak solution of L e {u e ) = div(/) + F mil. Then for any B C 2 B C fl, we have 
|Vti £ | £ L°°(B ) and t/ie uniform estimate 

||Vn £ || L oo( B ) < —|n e | 2 ) + C | ||/||l°o( 2 B) + r a [f] c o, a ^ 2B ) + |T P j , (4.3) 

where C depends only on [i,r, n,p,d,m and a. 

Proof. By rescaling we may assume r = 1, and £ (7“ (25) is a cut-off function such that (p = 1 on 5/4 B, 
ip = 0 outside 3/2 B, and |V</?| < C. Then we have 

-di v[AfV(<p U P)} = div(f a ip) + F a ip - f°Vip - AfVuPVip - di v(AfV^). 

It follows from the fundamental solution that for any x £ B, 

u £ (x) = - V y r £ (x,y)f(y)<p(y)dy+ T E (x,y)[Fip - f\7<p - A £ Vu £ Vp\dy +/ V :y T e (x, y)A E V<pu e dy, 

J2B JIB J2B 

(4.4) 

where we use integration by part in first and third term in right hand side of the above equality. Then 
differentiating both sides of (14.41) with respect to x gives 


Vu £ (x) = - / V x S7yT e {x,y)[f{y)ip{y) - f{x)tp{x)\dy - f(x)<p(x) / V x T £ (x,y)n(y)dS(y) 


I2B 


>d(2B) 


+ / V x T £ (x,y)[Fip-fVip- A £ Vu £ V(p\dy + / V x V y T £ (x,y)A £ V<pu £ dy, 


1 2B 


I2B 


where dS denotes the surface measure of dfl, and n is the outward unit normal to d{2B). We refer the 
reader to (22] pp.55] for the skill used above. 

Hence, in view of (|4.2[) . we obtain 


|Vu e (x)| < C< 

l 


\f(y)<p(y) - /(®M*) | 


I2B 


\x - y\ c 


dy + sup \f(x)ip(x)\ 

x£B 


\x—y\=l 1*^ V I 


dS(y ) 

d -1 


+ 


\ F (y)\ , 

k-y | d_1 


'(3/2B)\(5/4B) 


l/(y)l + |Vn e (y)| 
-id=i-^ + 


MiOl 


dy >, 


\ x _ y |d i 1 7 (3/2B )\ (5/4B) |® - y| d J 

where we use the observation of V<^ = 0 on 5/4 B and (p = 0 outside 3/2 B in last two terms. This leads to 
|Vrt e (x)| < C'|[/]c , o.<t( 2 _b) + ||/||l oo ( 2 B) + \\F\\ lp{2 b) + ( J |Vu £ | 2 dyj~ + j \u e \ 2 dy^~ 
for any x £ B. Then it follows from the Cacciopolli’s inequality that 


II Vu £ 11 loo( B ) < C{\\ u e\\L 2 {2B) + ll/llc°. CT (2B) + ll-f IIlp(2B)}j 

where p > d, and C depends on p,r, K,p,d,m and a. By using the rescaling technique, (14.31) can be easily 
derived. □ 
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Theorem 4.4. (Interior Lipschitz estimates). Suppose that A £ A (p,t,k), V satisfies (11.2)1 . (|1.4j) . and 
B,c satisfy (11.31) . Let p > d and a £ (0,1). Assume that u e £ F7/ oc (fl;M m ) is a weak solution to C £ (u £ ) = 
div(/) + F in LI, where f £ C 0,cr (f2; W nd ) and F £ L p (12;M m ). Then for any B C 2 B C Ll with 0 < r < 1, 
we have the uniform estimate 

||Vu e ||i<»(B) < — (^j |«e| 2 ^) _ + C | H/||l oo ( 2 _B) + r<7 [f]c°’ rT ( 2 B) + r (J~ i‘ ) | > (4-5) 

where C depends only on fi,r, k, X,p,d,m and a. 

Proof. We only need to prove (14.51) in the case of e < £o, where £o will be given later. Since the estimate (14.51) 
immediately follows from the classical results when e > eq. Consider the transformation T(x, e) = [T ^ 7 (x , e)] 
as follows 

= T /3l (x, e)vf = [8^ + ex q 1 (x/e)\v2- (4.6) 

In view of (12.131) . it is not hard to see T(x,e) is a diagonally dominant matrix whenever e < £\ = 
£i(y, r, k, m, d). Hence we have the existence of T” 1 (x, e), 

1/2 < ||T(-, e )|| ioo(n) < 3/2 and 2/3 < \\T~\; e)|| L<JO(n) <2 for £ £ (0, £l ), (4.7) 

where 

ll T (- ) £ )llL°°(n) = sup HT^j^lloo and \\T(x, e)||oo = max |T a/ 3 (x,e)|. 
xen l<a,p<m 

Moreover, in view of 


T-\x, £) - T-\y, e) = T~\x, e) [T(y, e) - T(x, e)] T~ 1 (y, e), 


we have 


|| T- 1 (x,£)-T- 1 (y,£)\\ 00 < \\T~\x,£)\U\T~\y,£)\\ 00 \\T(y,£)-T(x,£)\\ 00 < C\y - x\° 

for any a £ (0,1] and x,y £ B C fl, where we use [X'(-, e)]c7°>°'(si) — C'e 1_cr [xo]c 0 , tr (y) — C(h,t, n,m,d,a) 
which follows from (12.121) and (12.131) . Thus we obtain 

ll T_1 (-! £ )llco.‘Hfi) = l|F“ 1 (‘; e )IU°°(n) + [T _1 (-,£)] c0iCr(n) <max{2 ,C(y,r, n, a,m, d)} (4.8) 

for £ £ (0, £i). 

Consider the new system 

—div(A £ Vu £ ) = div(/) + F in f2, 

where Af 1 = Af 13 [5^ + £XqZ] > 

/« = /« + EV^Xolvf and F a = F Q + AfV X ^v] + V^VvJ - B^Vvg - (,cf + A 8 af3 )u 0 £ . 

Obviously, there exists £2 = £ 2 ( 44 , t, k. d, m) such that A £ A(L, r, k + 1) whenever £ < £ 2 . 

Let £0 = min{£i,£ 2 ,1} and £ £ (0,£o]. For any r £ (0,1], it follows from (14.31) that 

l|Vu e || L oo( S ) < cjr 1 1 |||£<x.( 2 ,b) + ||/||l°°(2B) + r<T [f]c°^'(2B) + r {J~ 2B l^l P ) P }’ (4-9) 

where a' = minjcr, r}. For convenience, we denote 

7 Z(nB) =-(£ |u £ ~) + ||/IU°°(n8) + r (^-[ |F/ J j J . 

J nB J nB 
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Hence, in view of (13.71) . (13.81) . (14.6j> and (14.81) . we obtain 
Ibelli 00 ( 2 S) < \\T 1 (■; e) ||l°°( 2 S) ll^e IIl°°( 2 B) < CrTZ(4B), 


\ v e]coy ( 2 B) — 1 (') e )]c°- CT '(2S) \\ u £\\l°°( 2 B) + 11^ 1 (’> e ) |Il°°(2B) [ u e}c 0 ^' ( 2 B) — C{r + F a }1Z(4B) 


- 1 , 


1 -a’ 


|Vu £ 


2 B 


< l|VT- 1 (-,e)|| L ^ ( 2 B) ||u £ || i ^ ( 2 i j ) + ||T” 1 (-,e)|| L ^( 2B )(^ |Vu £ |^x)* < C{r + 1}K(4B), 

(4.10) 

where we use Theorem EJ to estimate the term of (f 2B \X7u e \ p dx) l ^ P ■ 

Note that 

[f] C 0 ^'(2B) — (2B) + £ \yeXO,eVe\ C o (2B) 

< [/]c°’ CT '(2B) + e ll^llL 00 (R‘ 1 )IIXo||L°°(y)be]c 0 ^'(2S) 


+ £ 


l—o 


C 0, CT / ( K d)||Xo||l,°°(y)IKv||L°°(2B) + [Xo]cO, CT '(y)P IIl°°(R‘ 1 )II' (; £|Il°°(2B) 

— [f]c°' a ' (2B) “b ll 1 ’ellL 00 (2B) + [ v £}c°’ a '(2B) } 


} 


(4.11) 


where we use the condition (11.41) and the estimate (12.131) in the last inequality. Moreover, it follows from 
(14.101) and (14.111) that 

L°°( 2 B) + r a \f\c°’ a '{ 2 B) — ll/IU°°( 2 S) + r<7 [/]c°’ CT '( 2 S) + C\\ v e\\l°°( 2 B) + Cr a ||v e || c o,a'( 2B ) 

^ II/IIl°°( 2 S) + r<T lf]c°’ ,T '( 2 B) + C{r + r 1+a }TZ(4B) (4.12) 


and 


■(£ 


f i p r < 


'(£ 


< r 


2B 


2 B 


2 B 


F\ P Y+Cr(i |Vu £ | Py+CrU \Vu £ \ p ) P + Cr\\u £ \\ L 


2 B 


3 (2 B) 


F\ p ) p + C{r + r z }K{4B ) 


(4.13) 


Combining (14.91) , (14.101) , (14.12[) and ()4. 131) , we have 


II ||l°°(_b) < C’{ ||/||l°°( 2B) + rCT [/]c 0 > ct '(2B) + r(^j- l^l P ) P |+C'{l + r + r 1+<T + r"}7^(4H), 

where C depends on //, r, k, A, a,p, m, d. This, together with ( 12 . 121 ) . (14.71) and (14.101) . gives 

||Vu £ ||x,oo( S ) < ||VT(-,e)|| L oo( S )||u £ || I/ oo( S ) + ||T(-,e)|| i oo( B )||Vu £ || I ,oo( S ) < C{lZ(4B) + r a [f] c °,a( 4B )}. 


Note that since a 1 < a we have r a '[f] c — ^' ra [f]c°’ a (AB)- 

For any B r (x o) C B 2 r (xo ) C fi, there exist {Bl (xi)}^L 1 and x,; € B r (xo) such that B t {xq) C U F =1 Bi_(xi). 
It is clear to see B r (xi ) C H 2r (x o) for any x* S B t {xq). Hence we have 

l|Vu £ || L oo (SrOEo)) < max: {||Vu e || £ oo (Br(x . }) } < C{K(B r (xi)) + r a [f] c o,, {Br{xi)) } 


l<i<N 

c 


i r ~% 


1 

| 2 \ 2 


®2r(^o) 


+ 


>(B 2 r (xo)) + rCT [/]co. CT (S 2 r(a;o)) + r (/ I^T) P }> 

V B2r(x 0 ) ' J 


and we complete the proof. □ 

To prove the global Lipschitz estimates, we study some properties of the Dirichlet correctors 4> e fc, 
0 < k < d, which actually play a similar role as Xfc in the interior Lipschitz estimates. 
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Lemma 4.5. Suppose A £ A (g,T,n). Let g £ C' 0,1 (<9O; M m ) ; and u e £ i4 1 (r2;M m ) be the solution of 
L e (u e ) = 0 in Pi and u £ = g on dPl. Then for any Q £ dLl and e < r < diam(fi), 


' B(Q,r)nfl 


\\7u £ \ 2 dx 


<C{\\X7g\\ L °°(dfi)+ e 1 ||5 , IIl°°(9Q)} , 


(4.14) 


where C depends only on g, r, k, d, m and . 

Lemma 4.6. Let A £ A (g,T,n). Then we have 

\\®e,k ~ p kh°°(n) < Ce and ||V$ £ifc || L oo(n) < C (4.15) 


for 1 < k < d, where C depends only on g,r, K,d,m,g and fi. 


Remark 4.7. Lemmas 14.51 and 14.61 were proved in [38], as well as in [3], Here we omit the proof. 

Lemma 4.8. Assume that A £ A(g, r, k), and V satisfies (11.21) . (11.41) . Then for any a £ (0,1] , we have 

||$e,o — -f||L°°(n) < Ce, ||$ £ ,o — < Ce 1 a and ||V$ £) o||L°°(fi) < C, (4-16) 


where C depends only on g, r, k, d, m, g and Pi . 

Proof. Let ro = diam(H). If e > cro, (14.161) follows from the classical boundary Lipschitz estimates for 
elliptic system in divergence form with the Holder continuous coefficients. If 0 < e < cro, consider 


u £ (x ) = $ £ , 0 (x) ~ 1 ~ £ Xo{.x/e). 

Then L £ {u £ ) = L £ (4> £i o) — L £ [ex q(x/e)\ = 0 in fi, and u £ = —exo{x/e) on dPl. Hence, it follows from the 
Agmon-Miranda maximum principle (see [3j Theorem 3] or [38l Remark 3.4.4]) that 


sup \u £ (x)\ < C sup \u e (x)\ < C'e||xo||Lo°(Kd) < Ce, 
x€Q x£d£l 


which implies ||$ £ o — ^||L°°(n) A Ce. Additionally, for any a £ (0,1), in view of Theorem 11.21 we have 
ll* e , o — I\\c°^(Ci) A Ce 1 a . Note that L e is the special case of C £ , and C depends only on g, r, k, d, m, a and 

pi. 

Moreover, it follows from Lemma 14.51 that 


(-[ \Vu £ \ 2 dx) 1/2 < C, 

VD(Q,r) ' 


which implies 


(/ |V4> £)0 | 2 dxY /2 <C 

Vd(Q,v) ' 


for any Q £ dPl, and e < r < rg. By the interior Lipschtiz estimate, we have 

f 1 j 2 

sup |Vu e (a:)| < — (f \u e \ 2 dy) < C, 
{d x >e}nn £ V/o / 

which gives 


sup |V4* £]0 (x)| < C. 
{d x >e} nf2 

In the case of {d x < e} D H, we apply the blow-up argument. Let 


(4.17) 


(4.18) 


v (x) = ^ £;0 (ex) - 




then we have 
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Li(/u) = div(V) in fl £ , 

v = 0 on dLl e , 

where = {i £ : £i € 0}. Note that although the character of the boundary varies, fl £ is still a 

bounded C 1,v domain. The boundary functions of fl £ are denoted by tjji >e (x) = fifiex), i = 1, ■ • • , no (recall 
Remark 12.21) . and we fortunately have IIVvIIc 1 ^^- 1 ) £ e 1+r? ||?/ , j||c 1 -'' 7 (R d - 1 ) — £ 1+v Mq. Hence, it follows 
from the (boundary) Lipschitz estimate (|2.25l) that 

sup |Vv| <c\(-f |Vu| 2 dy) 2 + ||F|| c o, T(Kd) 

B(o,i)nn e I VB(o, 2 )nn s ' 



This implies 


sup |V$ £ , 0 (x)| <c\U |V4> £>0 | 2 dy) 5 + ||F|| c o,. (Kd) i < C. (4.19) 

x£D(0,e) ( VD(0,2e) 7 J 

Note that we choose r = 2e in (14.171) to give the last inequality. Combining (|4.18l) and (|4.19l) . we have 
l|V^ e ,o||z,<»(n) < C, and this also implies the second estimate of ()4.16D for a = 1. We thus complete the 
proof. □ 


Lemma 4.9. Assume the same conditions as in Lemma f478l Then we have 


||V4> e ,o|| c o,T(Q) < Cmaxje T ,l}. (4.20) 

Furthermore, exists and satisfies the following estimates: 

2/3 < s 2, ||V(#~2)|| i „ (n) < C, (4.21) 

whenever e < sq, where Eq = Eq (fj,,T,K,d,m,Ll) is sufficiently small, and C depends only on fi,r, K,d,m,r] 
and 14 

Proof. Let <h £j o = 4 > £i o — I, then L e (4> £i o) = div(H e ) in fl, and <l> £] o = 0 on <914 We first prove (14.20[) in the 
case of e < 1. Set U{e) =11(1 B(P, e) for any P £ 14 By translation we may assume P = 0. In view of the 
Schauder estimate (|2.26[) and Lemma 14.81 we have 

[V$ £i o] c o.r( [/ ( £ )) < Ce~ T j (/ |V$ £i o| 2 ^ 2 + ||14|| i oc(;7(2 £ )) +£ r [H e ]co,r( C 7(2 e )) 1 

( JU{2e) J 

< Cz r j||V<I? £i o||L°°(n) + l|l / llc°' T (R‘ 1 )} ^ Ce r , 

where C depends only on fi, r, k, m, d, rj and Mq. Note that [V$ £) o]co,i-(j 7 ( e \) = [V<1 ? £ ,o]c 0 - t (W e ))' Thus by a 
covering argument (see [22J pp.98]), we obtain ||V$ £! o||c 0 ' r (a) 4 Ce~ t . The case of e > 1 is trivial, since we 
can derive (14.201) by using the Schauder estimates (|2.19l) directly. 

Next we prove <|4.21[) . It follows from (14.161) that || < ho,e||L°°(o) < Ce. Since <l? £i o = I + 4 , o, e , we know that 
there exists 4>“o £ L°°(H) such that 

1/2 < || < 1? £ ,o||l°°(q) < 3/2 and 2/3 < ||^ £) ollL°°(f2) < 2 

whenever e < £o(p, r, k, d, m, H), and Eo is sufficiently small. 

Due to 

^o(*) - = ^o(*)[ $ e,o (y) - ^,o(^)]^o(y), 

we have 


|$ £ ,o(z) - $ e ,o(i0| < H^ e ,olU«»(n)||^ £j ollL«»(n)||V^ £ ,o||i,oo(n)|a: -y\<C\x-y 
for x,y £ LI, and this implies HV^olli/^Cn) — C*. The proof is complete. 


□ 
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Lemma 4.10. ( A nonuniform estimate). Suppose that A € A V satisfies (|1.2|) . (|1.4j) . and B,c 
satisfy (11.31) . Let p > d and a € (0, r]. Assume f € C' 0,o '(f2;M rnd ) and F € L p (fl;M m ) ; then the weak 
solution to C £ {u £ ) = div(/) + F in 12 and u £ = 0 on dPl satisfies the estimate 

l|Vu e || c o,f (n) < Cmax{ef- 1 ,l}{||/|| c o,a (a) + ||F|! LP(n )}, (4.22) 

where C depends only on [i,r, k, \,p,d,m,r] and PI. 

Proof. If e > 1, (14.221) follows directly from the Schauder estimate (|2.29l) and the Lipschitz estimate (12.281) . 
In the case of 0 < e < 1, the main idea is based upon the following interpolation inequality 

[Vfi e ]( q) < 2||V« e ||£ 00 ^ [Vu e ]^, 0iO .^. (4.23) 

Set U(e) = fill B(P,e) for any P € fl, and by translation we may assume P = 0. We first study 
l|V« e || L°°(n) through the uniform Holder estimates. To do so, let v £ = u £ — [I + exo{x / e))u £ {fi). Hence we 
have 


C £ {v £ ) = div(/) + F + [div(eI4xo,e) - B e V y xo,e ~ (c £ + A I)(I + £Xo,e)] «e(0) in PI, 
where y = x/e. If 0 £ dPi, we have v £ = 0 on dPi. Prom the Lipschitz estimate (12.311) on e scale, we obtain 


|Vu, 


'e ll-L 00 ((7(e)) ^ 


< 


< c{- 


1 

'e| 2 ) + II/I|l oo (C/( 2 £)) + e<T [/]c°' CT (t/( 2e)) + e 


U{ 2e) 


(/ |FP) P + M0)| 

v U (2e) 


m 


[fu 2 K -U e(°)| 2 ) 2 + C 'IKIU°°(fi) +C \\f\\c°^(n) + C£<7 \\ F \\lp(U( 2 )) 


’U(2e) 


< 


c 


'll Ue llc'°- <r (n) ^ {II ^ lichen) II-^IIlp(q)}’ 


- £ l-a ll u ' £ IIC°. CT (n) 
where we use ()3.14l) in the third inequality. This implies 

||Vn e || L oo(Q) < ||Vu e || L oo(Q) + ||Vxo||L° o (y)ll' u ellL°°(n) < Ce a 1 |||/||c ,0 ’ !T (n) + 11^11^(0)}• 
Next, it directly follows from the Schauder estimate ()2.32|) that 

1 


(4.24) 


[V 


U £\c°’°(U(e)) ~ C 


~1 ~\~(T 


1 

U e \ 2 ) + [f]c°> <r (U{ 2 e))+ e <T ||/IIl°°(C/( 2 £)) + £ 


U(2e) 


1—<7 


U( 2e) 


|FP 


< 


c 


— pl+cr 

c 


|||^e||.L°°(f2) + ll/llc'°.o-(0) + C*( j \F | 


1 

P\P 




where C depends on //, r, k, X,p, a,m, d,rj and II. By a covering argument (see [22j pp.98]), we have 

[Vu e ]co. CT (Q) < Ce 1 CT |||/||c 0 ' CT (a) + II^IIlp(Q)}- (4.25) 

Finally, the estimate (|4.22l) follows from (14.231) . (14.241) and (14.251) . We complete the proof. □ 

Proof of Theorem 11.31 In the case of g = 0, we only need to consider the following transformation 


«e(*) = ®e?o( x ) V 2( x 


(4.26) 
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for e < £*, where £* = min{£o,£ 1 ,£ 2)1 and £0 is given in Lemma H~9l and £i ,£2 can be chosen later. Since it 
is clear to see that the estimate (11.1011 immediately follows from the Lipschitz estimate (12.281) for £ > £*. 
Then the Dirichlet problem (11.71) can be transformed into 

f L e{v e ) = div(/) + F in 9, 

| v £ = 0 on 59, 

where 

f a = f a + Af{^ 0 -5^)Vv] + V e a ^^ 0 -8^)v2, 

F a = F a + + V £ a ^Vv] - BfX7u,P - (cf + X5 ap )u^. 

It follows from Theorem [14] and Corollary 13.81 that 

max{||Vu £ || L p (Q) , ||M £ || c o,a' (n) } < C{\\f\\ L oo m + H-FHi^n)} • 

This, together with Lemma 14.91 gives 

max{||Vu e || L p (Q) , |k||c°.-'(n)} < C {||/|| L oo (n) + ||^||z,p(Q)} , (4.28) 

where a' = 1 — d/p, and C depends on p,r, k, A, a,p, d,m and 9. Here we use v £ = $~q u £ , and Vu £ = 
V{<S>jl)u £ + d>”g Vu £ . However we need to rewrite Vu £ = —V<f> £i qv £ + (I — d > £j o)Vu £ + Vu E to handle the 
Holder norm of Vu £ . Set v = min{r, a, and v' = max{r, 1 — 2u}. Note that 0 < 1 / < 1, and we obtain 

[VUelcO.^Q) < [V^s^lcO^^llbelU 00 ^) + II V<l> ei o||L°°(fi) [Ve}c 0 ’-'(Q) + V ~ 4 > e ,o]c 0 ^(f2) II Vu £ ||x,oo(q) 

+ ||4 — 4> £i o||L°°(n)[Vu £ ]c'o, I ,(Q) + [Vu £ ] c o.,( Q) 

< C(e t + e 2v 1 ){||/||c'o,<t(q) + ||F||x,p(n)} + Ce 1 + Ce\^Jv £ \co^^ 

where we apply (14.161) . (|4.20l) . ()4.22 p and (|4.28p to the second inequality. This implies 

[Vu e ]( 70 ^(Q) < Ce v {H/llc 0 '' 7 ^) + II^IIlp(Q)} + C||Vu £ ||j / oo(q) (4.29) 

whenever e < £\, where E\ = min{l/(2C), 1}. Hence, we have 

WfWc 0 ’"^) < ll/llc°’ a (n) + C£||V« e |Uco (n) + [H e (d> £i o — -OVu £ ] co + IlK^^O — I) v e\\c°’ , '{n) 

< C'£ 1 -" , ||Vu £ || i oo (Q) + C{\\f\\ c o„ m + UTlI^p^)}, 

where we use (11.41) . (|4.16l) . (14.20p . (|4.28l) and (I4.29|) in the second inequality. In view of (I4.28p . we also have 

II-^IIl^o) < ll^lli,p(n) + C{ll^ e ||i,p(n) + ||Vu £ || LP (q) + ||u £ ||lp(q)} < C , {ll/llL oo (n) + ll-f Ilz,p(f2)}- 

We now apply Lemma 14.11 to (|4.27p and obtain 

IIVWeHicxqQ) < C’{||/|| c -o^(Q) + IlFUip^)} 

< C'£ W ||Vu £ || L 0 0(n) + C{||/|| c o, CT(n) + ||T IIlp(q)}, 

which gives ||V%|| L oo( n ) < C{\\f\\ c o,a (Q ) + ||F|| L p( n) }, whenever £ < £2 = min{l/(2C ), 1}. So we have 
||V« £ |U°°(n) < l|V( ( h £ ,ow £ )||L°°(n) < C{\\f\\ c o,<T + llTHip^)}, (4.30) 


where C depends only on p, r, k, A ,p, a, d, m, Mq,t] and |9|. 
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In the case of g / 0, consider the homogeneous system C £ {u £ ) = 0 in II and u e = g on 0, where 
g £ C 1,a (dfi;R m ) with a £ (0,?/]. Let h £ be the extension function of g, satisfying 

—div(A £ 'Vh £ ) = 0 in 0 and h £ = g on <90. 

It follows from Lemma ITT! that ||Vh e ||x,oo(Q) < C'||g , || f ;i, CT (,gQ), where C depends only on y,r, n,a,d,m,g and 
II. Let g = min{r, a}. By the argument applied to Lemma l4.10l we obtain [VL^o.^q) < Ce~ l ~ e ||5 , ||c 1 ’ CT (dr2)- 
Indeed, due to (12.2611 we have 

[ V M C o, W£)) < C£ ~ e {j B{p , ) l V ^l 2 Y ^ Ce-e\\Vh £ \\ Loom < Ce-°\\g\\ci,* iga) 

for any B(P, 2e) C O, while for the boundary estimates, it follows from the (boundary) Schauder estimate 
(10711 that 

[Vh £ \ C 0, e{ n {£)) — 7 ^{ll^ 7 ^ellL 00 (r>( 2 £)) + l|Vg||z / oo( A ( 2 £ )) +£ e [V5] C 0,« ? (A(2e)) +£ 1 IIS , llL°°(A( 2 e)) j 

< ^l^{ll V ^elk°°(n) + IMIc^^n) j < Ce 1 e ||5'llc 1 ’ CT (90) 

for any P £ <90, where C depends on y, r, k, q, d, m, Mq, rj, and 101. Thus we have 

I i i+„ 

[Vh £ ] c o, f(c) < 2||Vh B ||>. (n) [V^]> 0i<l(n) < Ce-—\\g\\c^ {d ny (4.31) 

Set We(x) = u £ (x) — $^g(i)/ie(x), we obtain 

( £e(w e ) = div(Z) + F, in 0, 

w £ = 0, on <90, 

where 

f a = Afi^l-d^VhJ + V^i^-d^h], 

F a = AfV^Vh] + V^VKl - BfV^lKl) - (,cf + Xd^^h]. 

Now, let v = g/2. In view of (|4.16H and (14. 3111 . we have 

ll/llc°^(n) < C'll5 , llc 1 - CT (9n) and H-^llLP(f 2 ) < C'II5 , IIc 1 - ct (9o)- 

Note that ||||i°°(f 2 ) < C\\y\\L°°(dCi)i where C depends only on y, r, k, er, d, rn and 14. Hence, recalling (14.30(1 . 
we have 

||Vtl e || L oo(Q) < HVlUellioo^) + j|V(<h £i oh £ )|| L oo(Q) 

< C{\\f\\c°^(n) + ll^llLP(n) + lls , llc 1 - CT (dQ)} ^ C\\g\\ci.*(dn)- (4.32) 

Finally, (11,101) follows from (14.3011 and (14.32(1 by writing u £ = u e> i + 2 , where u £ y,u E} 2 respectively 

satisfy the homogeneous and non-homogeneous systems (see (13. 151) 1. The proof is complete. □ 

Lemma 4.11. Suppose that the coefficients of C £ satisfy the same conditions as in Theorem, II. 41 Then 
G e (x,y) has the following estimates: 

\VxV y Ge{x,y)\ < ^ _^| d , (4.33) 

\VxGe{x,y)\ < j^-~pr m in{l, and I VyGe(x,y)\ < minjl, ^ } (4.34) 

for any x,y £ Q and x y, where C depends only on y, r, n, A, d, m, g and 0. 
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Proof. For any x,y £ O, let r = \x — y\. Due to 

£*e[Ge(-,y)\=0 in fl\B(y,p) 
for any p > 0, it follows from (11.1011 and (13.321) that 

C 


^ vflr ( X ) 


\VxGe(x,y)\ < ||Va e (-,2/)|| L oo (f 2r(a:)) < — ( i I G £ {z,y)\ 2 dz\ ‘ < Cr 1 d , 


(4.35) 


(4.36) 


where x can be on <90. By applying the localization technique (as shown in Remark l2.13l) to (|1.10l) . we have 

1/2 


||Vu £ ||L°°(nr(i)) < — (£ 


(x) 


2 dy ) 


for u £ satisfying C e {u £ ) = 0 in On (a;) and u £ = 0 on <9(Oi(x)) D <90. (We remark that we just consider 


the estimate at boundary, and the interior one directly follows from (14.311 .1 So, we can derive the second 
inequality of (|4.36l) . 

For the adjoint Green’s matrix *Q e (-,x), we have \\7 y Q £ (x,y)\ = |V*C/ e (y, x)| < Cr 1 ~ d by the same 
argument. Moreover, since S7 y G e (-, y) still satisfies (j4.35f) for any p > 0, and S7 y G £ {-, y) = 0 on <90, we obtain 


\VxV y G e (x,y)\ < 


C 


(£ 


fir (#) 


\z - y\ 2{l ~ d) dz ) 2 < Cr 


n —d 


where r/2 <\z — y\< 2r. Observe that V y Q e (-,y) = 0 and S7 x [*G e {’,x)] = 0 on <90, we have 

NyG £ {x,y)\ = |Vj ,G e {x,y) - V y G £ (x,y)\ < \V x V y G e (x, y)\\x - x\ < 

where x £ <90 such that d x = \x — x\. Similarly, we have \S7 x G e {x,y)\ < Cd y \x — y\~ d , and the proof is 
complete. □ 

Proof of Theorem 11.41 Define the conormal derivative (J^)* corresponding to C E and C* as 
follows: 


d d 

— = -niA i:j {x/e)— - mViix/e), 


( d d 

( QGj = -^^(x/e)— - njBjix/e), 


where n = (ni,--- , n^) denotes the outward unit normal vector to <90. Thus, define the Poisson kernel 
V £ (-,y) = [V] ,3 (-,y)\ associated with C £ as 


n\x, y ) = 


d_ 

dv 


[G £ ^(x,y)] = -n j (x)a?f(x/e)-^-{g? r (x,y)} -n j (x)B* fi (x/e)gf 1 {x,y) 


for y £ 0 and x £ dfl. It follows from (14.3411 that 

\V £ (x,y)\ < Cd y \x - y \~ d , 


(4.37) 


where C depends only on p, r, k, A, d, m, rj and O. Thus for any g £ L p (dfl; M m ) with p £ (1, oo], the solution 
to C £ {u £ ) = 0 in O and u £ = g on dfl can be written by 


u s (y)= V £ (x,y)g(x)dS(x) 
J an 


for any y £ fl, and it follows from (I4.37p that 

\u e (y )| < Cd y 


Ian 


\g(x)\ 

\x - y\ c 


dS(x) 


(4.38) 


(4.39) 
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Recall that the nontangential maximal function of u £ is defined by 

(■ u £ )*(Q ) = sup{|u £ (x)| : x £ n and \x — Q\ < N 0 dist(x, (90)} for Q £ dLl, (4.40) 

where Nq = A r o(D) > 1 is sufficiently large. Hence, if \y — xo| < N$d y for some xq £ dfl, then we have 

\u e (y)\<cl \g(x)\dS(x) +Cdy^2 [ 

J dttnB(xo,r) j—QJYij V I 


oo p 

<c^+ \g(x)\dS(x ) + Y] 2 \g(x)\dS(x] 

Jdi 


/ 

J dflnB(xo,r) 


< C sup 

0<r<diam(f2) dfi,r\B(xo,r) 


■ = q o dQrB(x 0 ,2o+ 1 r) 

\g(x)\dS(x ) 


where r = d y /2 and £j = 90 n {l?(xo, 2- J+1 r) \ B(x o, 2 J r)}. Note that 

•A^9n(bl)(^o) = sup {/ 

0<r<diam(Q) d£lC\B(xo,r) ^ 

is the Hardy-Littlewood maximal function of g on <90. Thus it is not hard to see that 

(u e )*(x 0 ) < CM d n{\g\){xo). 

Due to the L p bounded properties of the Hardy-Littlewood maximal operator: ||-A^ar 2 (|<?|) Ilz,p(ar 2 ) < 
C\\g\\LP(dn) for 1 < p < oo (see [41] pp.5]), the estimates (I1.14J) and (11.151) can be derived immediately. 

Now, we turn back to verify (I4.38p . Let R = d y / 2, and <p £ Cq{B(jj,R)) be a cut-off function such that 
(p = 1 in B(y,R/ 4) and ip = 0 outside B(y,R/ 2). Then, since C £ {u £ ) = 0, we have C £ (<pu £ ) = —C £ [(l — ip)u £ \ 
in O and pu £ = 0 on <90. Hence, in view of (I3.18p . we obtain 

(pu £ )(y) = - [ g £ (-,y)C £ {{\-y)u £ \ 

Jn 

= l Ge(-,y)£e[(l - <p)u £ ] - [ C* e [g £ {-,y)\[(l-cp)u £ }\ 

r ^° 1 Jn\B(y,r) Jn\B(y,r) J 

= I 'Pe(-,y)[( 1 -v)ue] + lim [ - p)u £ \ - lim [ V £ (-,y)[{l - <p)u £ ] 

JdU r->-0 JdB(v.r) 6lJ r ^°JdB(v,r) 


[ 

Jan 


'Pe{-,y) K 1 - <p)u e \. 


Note that C*[Q e (-,y )] = 0 in fi \ B(y,r ) for any r > 0, and (1 — <p)u £ = 0 in B(y 1 R/ 4). The proof is 
complete. □ 

Remark 4.12. Note that the same type of results for L e with Dirichlet boundary conditions and Neumann 
boundary conditions were shown in [5] Theorem 3] and [311 Theorem 1.3], respectively. Also, we refer the 
reader to [2j Theorem 1.3] for the same type of result in the almost periodic setting. In the case of m = 1, 
when we derive the estimate (11.151) with (7 = 1, there is no regularity condition on the coefficients of C £ , 
but some additional conditions on V are inevitably required even when A > Ao, and e = 1 (see [22] pp.179]). 

5 Convergence rates 

Lemma 5.1. Suppose that A £ A(h,t,k), and V satisfies (11.21) and (11.41) . Let 

= $“o( x ) - ~ TAoV/e), t( x ) - P k - £ xl{.x/e). 
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Then we have 


max ||^ £ / c || £ , t x>(Q) < Ce, max |V^ £ifc (x)| < Cminjl, ed^. 1 }, 

0 <k<d o <k<d 


(5.1) 


where C depends only on fi,r, K,m,d and O. 

Proof. By the definition of Dirichlet correctors & £>k , we have L e {}.V £tk ) = 0 in 11 and = ~ £ Xk,e on dO. 
Thus it follows from the interior Lipschitz estimate ()4.3I) and Agmon-Miranda maximum principle (see [3) ) 
that, 

|VT £j fc(x)| < \^ e ,k\ 2 dy) 1 < Cedf 1 , Vx £ 0. 

The rest parts of the lemma follow from Lemmas 14.61 and 14.81 and Remark 12.91 □ 

Lemma 5.2. Suppose that u e £ iL 1 (H) J u £ H 2 (0 ) and C s (u e ) = Cq{u) in 0. Let 

, du 1 


We =U P £ ~ $f,> 7 - j. - X k 5^] 


dxP 


where 1 < k < d. Then 


d 


[£e(u> e )] {^ k ^dx h dx 


+ F, 


i cry 


du 7 


kOXj 

K j 


kl,£ dx k 
d 2 u< 


+ HT1uA + 


a 0 \ d *el 9 2 VT t 9< 7 0 du^ 


hi 


l i ,£ ] dxj dx k dxi 


+ 


dxidx k 


+ [<5 - «*] 


du 7 
dxi 


3*1 J 1 

*’ £ | dxi dx k dxi J 
dd?F d 2 rf dC, ai dtf 


dxj dxi 
hi 


- B 


[cf + \ [®% - Xt s*] fL + _ 5 fh] u 


7 1 - 


— £- 


dyt dx k dxi dyi dx 


(5-2) 


where y = x/e, and 


_ 


^ + <&[*R-* k * fh ] 


ni = **£ + < £ [< 7 - <^ 7 ] + [< 7 - x fc d^] + e . 


ddz 1 

dyi 


Ke 


07 = c«/[ 4 .fj-^]+ £ A!!. 

fVote that ,, F(( 7 , d)T', are defined in Remark [2791 
Proof. From C e {u £ ) = £o(^)> it follows that 


clu 


[4(%)]“ = [£o(«)] a - [-C £ (4 > £ ,o«)] q - [£e((* e ,fc - x fc /)—)] 


{^Sj} + + P" + A«“V - - N“. (5.3) 


_d_ 

dxi I cte 


where h = £ £ (<F £i0 w) and J 2 = £ £ [($£,* - x k I)^\. By the definition of < f> £i fc, and Xki 0 < k < d, we obtain 


[h] a = 




duP 


d 


du[\_ a0 g < 7 o du'y 

dx,, / ai ^ £ dxn dx; 


8 Ks} [< 7 o - + Bg A.{<>3} + [ c f + A 


dx 
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d 


a(3 


duP 


\- 


d 


dxi l u,£ dxj J dx. 


I af lx.Pl sPl\ 8 “ T 1 of Bill „,j flXo 1 8u' 


sft? 1*3 * **]« 7 } - Kf ^+ [*5 * ^ ■ + ^“ 7 + 1H 


((Tu^ 


(9u^ 


<9u 7 


d Vj 

3Tf 7 


<9x; 


dy 


fii 

^u 7 


+ [cf + A<5 Q/3 ]$f 7 u 7 


and 

fc] a 


dxil ij ’ £ d Vj dxj 5x,;l ^ £ <9x,- L £ ’ fc fc Xfe - £j «9x fc J 5x,;l^’ £ ^ £ ’ fc k ^ dxj 


dxi 

d 2 ui 


U<^ - ^ S + <U [*S - ^ £}+i * + >**\ [*S - ***] £ 


<9x fc . 
9x fc . 


Sxfc 

du 1 


d f a/3 d Xj 7 d^ 7 ) _ a/3 

8iil %f 9® 3iJ 


8*fJ 


<9 2 « 7 


3 


<9x 7 - dxkdxi dxi 




<9 2 ii 7 

dxjdxk 


} 


_d_ 

dxi 


l * ,£ ^ £,fc ^ J l,£ l dyi dxk dxi dxi, £,fc dxjdxk J 


<9u 7 


+ [c^ + Ar^Kj-x^]^, 

where y = x/e. Put ii and I 2 into (15.31) . and then we have 


<9 2 u 7 

dxkdxj 


+“^[<5-^1^} 




+ ft, 


<9u 7 


_ B afi I - x,^ 7 i d2ul + - sM — + 

*’ £ I ™ k J 3®^®*. + ^ £ ’° J dxi + dxi 8x k + dxi 


-,M 


- [<P + 


<5 - + [<I - ^if?} - + zr<(s,4) 


<9u 7 


,<9u 7 


J 


< 9 ®,; 


dan 


where 6“ 7 , £^ a7 , VP. 0 ) 7 , are defined in Remark 12.91 Besides, the following identities hold. 


-£K<} - -£{£raj£} 




d 2 u J 
£ dx k dxj 


= e-£- <! p;“2 

axi 


uTiv) 


w, 


aj 


du 7 
<9xj 
9u 7 






0Xj 


dx k 


ai du 7 X 
fc * ,£ dxi J 


9 2 u 7 

kl3 ' £ dxkdxj 
du 1 


cP*Vb^ 3 f 

— eF ai = c < F a7 _ 

fc * ,£ dxkdxi dxk \ kl,£ dxi 


k ’ £ dx k 
Z^u 1 


c _d_ (dtf^duZ 

dxi \ dyi dx k 


= £ 


— £ 


dft^i q2 u j 


d f d( ai ^ 

£ dxi { d Vi U £ ' 


dyi dx k dxi ’ 

<9C“ 7 du 7 


dxi dxi 

These together with (15.41) give the formula (15.21) . and we complete the proof. 

Proof of Theorem 11.51 Let w E ^,w E ^ £ iLo(fl;M m ) satisfy w E = + w^ 2 , such that 


□ 


£e{We,l) = £e{We) ~ © in fl, 
where w £ is given in Lemma 15.21 and 0 = (0“) satisfies 


C £ (w E o) = © in fi, 


(5.5) 
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0" = a afi 


d^el d 2 v? d^'ll dvB 

dxj dxi 


lJ ' £ I dxj dxkdx 


- B 


0.(3 


' du< 


l,£ I dxi dxk 


+ 


dxi U 


For the first equation of (15.5p , it immediately follows from Theorem ll.il Lemmas 14.6114.81 and Remark 12.91 
that 


l|Vw ej i||ip(f2) < Cs\\u\\ W 2,p^y (5.6) 

For the second equation, in view of (12.71) . we have 

Col|Vu; £> 2|li2 (n) < B £ [w £t2 ,W £j2 ] = J Q< * W £ 2- 

For the right hand side, it follows from (|5.1I) and Cauchy’s inequality that 

[ ©“<,2 <Ce [ (|V 2 u| + |Vu| + \u\)\w £ , 2 \ ^ < Ce\\u\\ H 2 (n) \\Vw £)2 \\ L 2 (sl) , 

Jn Jn ix 

where we use Hardy’s inequality in the last inequality. Hence we have 

II^H^CO) < ^{ll^^e.illiRn) + l|Viu £ ,2|| L2( n)} < Cs\\u\\ H 2^, 

and this gives the estimate (11.161) . 

We now turn to show the estimate (|1.15l) . First, by recalling the estimate (|3.33l) . we have 


\Oe(x,y)\ < 


Cd y 

x — y\ d ~ l 


(5.7) 


Additionally, in view of (15.11) . we have 


16(9)1 < T-{|V 2 a(9)| + |V»fe)| + |u(„)|} (5.8) 

for any y £ H. 

Since 

w £ 2 {x)= [ g £ (x,y) a PQ a (y)dy for x € fi, 

Jn 

it follows from (15.71) and (15.81) that 

\we, 2 ( x )\ < Ce J Q , | a; _y|d-i '{l v2tt (y)l+ l Vu (y)l + \ u (v)\} d y- (5.9) 

Thus by the Hardy-Littlewood-Sobolev theorem of fractional integration (see m pp.119]), we obtain 


< CE\\u\\ W 2, P (n) 

with 1/q = 1/p — 1/d when 1 < p < d. For p > d, we can straightforward use Holder’s inequality to arrive 
at 


||w£,2||L°°(fi) < Ce\\u\\ W 2, P ^y 

Besides, it follows from (15.61) and the Sobolev inequality that ||io e) i||L«(fi) < C|| ||z,p(o) < C£\\u\\w 2 ’P(n)- 

Thus we obtain 

II^eIIlRQ) < {||^e,l||L«(n) + II’^e, 2 ||l 9(fi) } < C£\\u\\ W 2 ,p(ty, 
which implies the estimate (11.171) . and the proof is complete. □ 
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Remark 5.3. In view of Lemma 14.111 and Lemma 1 5.2 1 we can actually derive the following estimates by 
the arguments developed in m , 


G e {x,y) - Go(x, y) 


Ce 

x — y| d_1 ’ 


V x,y € n and x / y. 


Then we have 

IK -«IIl°°(Q) < Cefln^oe -1 + 2)] i ^ 3 ||F|| Z/d(0) , 

where Rq denotes the diameter of fL Moreover, let fl be a bounded C 2,ri domain, we have 


< 


Celn(e l \x — y \ + 2) 
\x — y\ d 


for any x, y € 11 and x ^ y. Then it follows that for any 1 < p < oo, 

IK - *./»- (<* - P?)^- k \\wi-m S CellnCRoe- 1 + 2)] 2| K l ||U|„ (! r l , 

where C depends only on //, r, k, A, m, d,p and fl. The details are left to readers (or see 
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